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EXECUTIVE  SUMMARY 

A critical  part  of  the  analysis  and  estimation 
of  sonar  system  performance  is  the  characterization  of 
the  acoustic  signal  and  background  noise.  Current  tech- 
niques in  performance  and  engagement  modeling  employ  random- 
processes  to  simulate  the  long  and  short-term  variations 
of  these  parameters  in  time.  The  usual  justification  is 
that  the  complex  mechanisms  governing  the  generation,  trans- 
mission, reception,  and  processor  response  of  signals  and 
noise  are  sufficiently  "random"  to  be  modeled  by  stochas- 
tic processes  determined  by  a few  input  parameters.  Since 
such  an  approach  is  much  more  efficient  than  measuring  or 
modeling  the  details  of  the  acoustic  environment,  it  is 
widely  utilized.  However,  it  has  not  undergone  a compre- 
hensive evaluation. 

This  report  documents  an  investigation  into  sev- 
eral aspects  of  the  validity  of  the  random-process  approach. 
Among  the  topics  considered  are: 

(a)  the  basic  statistical  properties  of  acoustic 
signal  and  noise  fluctuations 

(b)  the  ability  of  the  random-process  models 
to  properly  simulate  the  important  proper- 
ties of  the  acoustic  variables 

(c)  the  ability  of  the  stochastic  models  to 
accurately  predict  measures  of  system  effec- 
tiveness 
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(d)  the  choice  of  a stochastic  model  and  esti- 

! 

mation  of  its  inputs,  for  a particular 
scenario . 

The  ideal  way  to  address  these  would  be  to  compare  stochas- 
tic simulations  with  measured  acoustic  and  performance 
data  for  many  ocean  environments,  sonar  systems,  and  en- 
gagement geometries.  However,  the  scope  of  this  study 
and  the  availability  of  appropriate  measurement  data  re- 
quired the  investigation  to  focus  on  a simplified  test 
case  for  a representative  acoustic  environment,  sonar 
system,  and  target. 

A single  ocean  environment  (in  the  Northeast  Pacific 
dviring  Summer)  and  target  scenario  has  been  selected  for 
the  test  case.  To  ensure  that  the  effects  of  the  environ- 
mental acoustic  parameters  could  be  gauged,  the  study 
assumed  a stable  narrowband  source  (at  25  Hz),  an  idealized 
towed  array  receiver,  and  a processing  system  which  auto- 
matically detects  when  the  beam-output  signal-to-noise 
ratio  (SNR)  exceeds  a threshold.  Evaluation  of  the  random- 
process  approach  under  these  conditions  then  requires  the 
detailed  space/time  properties  of  the  transmission  loss 
(TL)  and  ambient  noise  A novel  aspect  of  the  present 
investigation  is  the  use  of  state-of-the-art  acoustic 
propagation  and  ambient-noise  models  to  fabricate  such  a 
controlled,  representative  acoustic  environment. 

For  the  test  scenario  and  target  conditions, 
there  is  good  reason  to  attribute  fluctuations  in  TL  to 
relative  receiver  movement  through  the  source's  multipath 
interference  field  (assuming  little  influence  from  ocean 


dynamics).  The  acoustic  TL  "data"  were  generated  from  mo- 
dels which  fully  account  for  this  mechanism.  Similarly, 
a special  ambient  noise  model  was  developed  to  simulate 
the  dominant  noise  fluctuation  mechanism:  movement  of 
many  contributing  surface  ships  through  the  interference 
field. 

Once  signal  and  noise  time  series  were  constructed 
with  the  acoustic  models,  several  ot  the  more  popular  ran- 
dom-process models  were  used  to  simulate  the  series.  In- 
put parameters  (the  trend  in  TL,  the  variance  and  decorre- 
lation scales  of  TL  and  noise)  were  determined  directly 
from  the  acoustic  data.  Comparisons  were  then  made  of 
the  statistical  properties  of  the  acoustically  and  stochas- 
tically modeled  signal  and  noise,  as  well  as  of  the  detec- 
tion properties  of  the  signal-to-noise  ratio. 

The  principal  conclusion  of  the  study,  subject 
to  the  limitations  in  the  test  environment  and  data,  is 
stated  generally  as; 

e Given  proper  input  parameters,  the  stochastic- 
process  approach  can  produce  adequate  simu- 
lations of  signal  and  noise.  The  accuracy 
is  sufficient  to  provide  useful  estimates 
of  detection  probabilities  and  detection 
histories. 


This  result  is  tempered  by  the  observations  that: 

• The  accuracy  and  utility  of  the  approach 
are  greatly  dependent  on  the  availability 
of  accurate  input  data  (e.g.,  the  variance 
of  beam-noise  fluctuations).  Such  data 
may  be  difficult  to  obtain,  and  without 
them  the  predictions  will  usually  be  poor. 

• The  accuracy  also  depends  on  the  particular 
random  process  used.  This  study  concentrated 
on  the  popular  ones,  and  the  agreement  with 
acoustic  data  was  just  adequate.  However, 
significant  improvements  can  be  realized 
when  the  basic  statistical  properties  (e.g., 
marginal  density)  are  tailored  to  the  vari- 
able and  case  at  band. 

The  report  presents  the  details  of  the  comparisons  on 
which  these  Judgements  are  based,  and  discusses  the  esti- 
mation of  input  data,  the  choice  of  a particular  random- 
process  model,  and  the  prediction  of  measures  of  effective- 
ness . 


Finally,  recommendations  suggested  by  the  inves- 
tigation include 

• use  of  the  computer  packages  and  Insight 

developed  for  the  study  to  analyze  measure- 
ment data  and  to  extrapolate  the  results 
to  other  environments,  receivers,  and  en- 
gagement scenarios. 


• application  of  the  theory  and  analytical 
models  for  multipath-induced  fluctuations 
and  ambient  noise  to  generate  better  random- 
process  models  and  inputs. 

The  details  of  the  study  are  presented  in  a two- 
volume  report.  This  first  volume  summarizes  the  modeling 
procedure  and  gives  the  results  and  conclusions  of  the  com- 
parison of  the  stochastic  process  results  with  the  acoustic 
data.  The  companion  volume  is  a compilation  of  appendices 
which  provide  references,  calculations,  documentation  and 
background  information  to  support  the  results  reported  in 
Volume  I . 
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Section  1 

INTRODUCTION  AND  SUMMARY 


1.1  The  Problem 
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Present  techniques  for  predicting  and  evaluating 
the  performance  of  Navy  sonar  systems  depend  heavily  upon 
stochastic  models  to  generate  basic  time  series  of  signal 
and  noise  or  signal-to-noise  ratio.  In  particular,  the 
deviation  from  a deterministic  mean  value  is  most  often 
treated  as  a random  process  with  specific  temporal  (e.g., 
stationary,  Markovian  with  exponential  autocorrelation 
function)  and  distributional  (e.g.,  Gaussian,  Rayleigh, 
Rician)  properties.  Time  series  of  signal-to-noise 
ratio  are  thus  generated  as  realizations  of  the  random 
process  and  passed  through  detectors  to  simulate  detec- 
tion histories  for  individual  encounters.  Statistical 
accvimulation  over  a number  of  such  encounters  yields  esti- 
mates of  probabilities  of  detection  and  other  measures 
of  system  effectiveness,  such  as  "mean  holding  time" 
and  "time  to  lost  contact." 

The  critical  assumption  in  such  an  approach  is 
that  those  aspects  of  the  slgnal-to-noise  environment 
which  control  the  detection  process  are  being  accurately 
simulated  by  the  specific  stochastic  process  with  its 
assumed  parameters.  The  usual  Justification  is  that 
the  extremely  complicated  processes  governing  the  genera- 
tion, transmission,  and  reception  of  acoustic  signals 
and  noise  in  the  ocean  are  sufficiently  "random"  to  be 
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■ characterized  in  terms  oi  a stochastic  process  which  de- 
pends on  a limited  number  of  parameters  and  is  for  the 

? most  part  independent  of  the  details  of  the  acoustic  en- 

■ ^ vironment.  Clearly,  if  this  is  the  case,  such  simula- 

tion techniques  offer  a considerable  savings  in  effort 
over  the  measurement  or  modeling  of  the  entire  acoustic 
environment  in  terms  of  the  underlying  transmission  loss, 

% noise,  target  signal,  sonar  response,  etc. 

f This  study  then  addresses  the  question  of  val- 

idity of  the  stochastic-process  approach;  its  ability 
to  simulate  the  important  properties  of  the  acoustic  en- 
vironment and  to  provide  accurate  estimates  of  system 
performance.  Implicit  in  this  question  are  such  problems 
as  the  choice  of  a particular  random-process  model  for 
a given  scenario,  the  determination  of  its  input  parame- 
ters, and  the  selection  of  a technique  to  actually  gener- 
ate sample  time  series  of  signal-to-noise  ratio. 


Ideally,  the  validity  of  the  random-process 
approach  might  be  determined  for  a particular  type  of  sonar 
system  and  engagement  as  follows . 


(a)  Select  a set  of  environments  which  spans 
conditions  of  interest  for  the  type  of 
engagement  and  sonar. 

(b)  Carefully  execute  at-sea  experiments  to 
measure  the  detailed  signal  and  noise  pro- 
perties for  each  environment,  and  for  a 
variety  of  source/receiver  geometries 

and  movements.  Use  a controlled  source 
and  the  selected  sonar  system  as  receiver. 
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(c)  Subject  the  beamformer  outputs  to  the 
sonar's  processing  and  detection  systems. 
The  collected  data  would  then  consist  of 
time  series  of  signal  and  noise  and  de- 
tection for  a variety  of  environmental 
and  source/receiver  conditions. 

(d)  Select  random-process  models  and  asso- 
ciated input  parameters  to  simulate  fluc- 
tuation and  detection  properties  for  the 
signal,  noise,  and  detector  input  (e.g., 
signal -plus-noise) . 

(e)  Generate  statistical  properties  of  sig- 
nal and  noise  with  the  random-process 
models  for  the  conditions  of  the  experi- 
ments. Compare  these  simulation  results 
with  the  measured  data. 

(f)  Subject  the  random-process  signal  and  noise 
data  to  the  sonar's  processing  and  detec- 
tor algorithms.  Compare  the  resulting  de- 
tection histories  and  other  measures  of  ef- 
fectiveness with  those  observed  in  the  mea- 
surements. 

(g)  Judge  "validity"  of  the  random-process  ap- 
proach from  the  accuracy  of  the  simulated 
signal  and  noise  properties  and  detection 
histories. 


Each  of  these  steps  is  certainly  feasible,  but  the  experi- 
mental effort  would  be  overwhelming,  even  if  only  a single 
environment  were  treated. 


In  the  present  study,  the  seven  steps  are  fol- 
lowed, but  with  alternatives  for  (a),  (b),  and  (c)  which 
allow  the  problem  to  be  attacked  with  limited  resources: 


(a)  Perform  the  analsls  for  a single  environ- 
mental case. 

(b)  In  place  of  at-sea  measurements,  construct 
representative  acoustic  data  with  acoustic 
models.  Limit  the  number  of  source/recei- 
ver geometries . 
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(c)  Use  an  idealized  sonar  system,  with  time- 
invariant,  deterministic  beamforming,  pro- 
cessing, and  detector  algorithms. 

The  choice  of  (a)  limits  the  ensemble  of  conditions  treated 
here,  a constraint  imposed  by  the  scope  of  the  study.  Use 
of  an  idealized  sonar  in  (c)  not  only  simplifies  the  data 
analysis,  but  also  allows  the  study  to  concentrate  on  the 
environmental  acoustic  fluctuations  unperturbed  by  uncer- 
tainties or  variations  in  sonar-system  behavior. 

The  most  important  substitution  (b)  is  of  acous- 
tically modeled  data  for  measurement  data.  Cost  savings 
alone  make  it  attractive.  But  it  also  provides  the  oppor- 
tunity to  investigate  the  sensitivity  of  the  fluctuations 
to  the  underlying  environmental  and  geometric  factors  by 
isolating  each  influence  and  identifying  the  dominant 
mechanisms.  Confidence  that  the  modeled  data  are  realis- 
tic follows  from  recent  advances  in  low-frequency  propa- 
gation and  noise  modeling  in  which  it  has  been  suggested 
that  not  only  averaged  transmission  loss  (TL)  and  noise, 
but  also  their  fluctuation  properties  can  be  accurately 
predicted.  Modelable  characteristics  include  overall  dis- 
tributions of  intensity,  range  spectra  and  arrival  angles 
for  TL,  time  spectra  and  directionality  for  noise.  This 
capability  makes  it  feasible  to  create  representative 
signal  and  noise  environments  with  realistic  fluctuation 
properties . 

In  summary,  this  study  concentrated  on  a single 
North  Pacific  oceanographic  environment,  a few  source 
tracks  relative  to  a fixed  receiver,  a stable  narrowband 


t 


1-4 


source,  a generic  type  of  towed  array  system  and  an  auto- 
matic detector  operating  on  signal-to-noise  ratio  (SNR). 
Signal  and  noise  time  series  were  generated  with  state- 
of-the-art  acoustic  models  for  the  idealized  sonar  re- 
sponse function  and  target.  The  acoustic  model  for  trans- 
mission includes  what  is  believed  to  be  the  important 
fluctuation  mechanism  for  the  test  case  (multipath  inter- 
ference for  a moving  source,  but  for  a static  medium). 

The  automatic  detector  and  statistical  analysis  routines 
were  applied  to  the  signal,  noise,  and  SNR  series  to 
yield  statistical  properties  of  the  acoustic  variables 
and  detection  histories  for  the  various  target  tracks 
and  sonar  characteristics.  These  are  the  "acoustically 
modeled"  data. 

Given  the  acoustic  data,  selected  statistics 
of  the  signal  and  noise  time  series  served  as  input  to 
stochastic  models,  which  simulate  fluctuations  by  treat- 
ing them  as  parameterized  random  processes . The  focus 
was  on  the  Gauss-Markov,  Jump,  and  Ehrenfest  processes, 
each  of  current  interest  in  Navy  operations  analysis 
and  determined  by  at  most  three  parameters.  "Stochasti- 
cally modeled"  data  were  thus  generated  for  the  test 
case . 

Comparisons  of  the  acoustic  and  stochastic 
data  were  then  made  with  the  following  questions  in  mind: 

(i)  What  are  the  basic  statistical  properties 
of  acoustically  modeled  signal  and  noise 
fluctuations? 
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(ii)  Under  what  conditions  can  the  stochastic 
simulations  properly  model  the  signal  and 
noise  time  series? 

(iii)  Under  what  conditions  can  the  stochastic- 
process  models  accurately  predict  typical 
measures  of  system  effectiveness? 

(iv)  How  can  predictive  acoustic  modeling  be 
used  to  aid  in  the  selection  of  random- 
process  models  and  parameters  for  specific 
engagement  scenarios? 


1 . 3 Summary  and  Outline  of  Report 

The  results  of  this  investigation  are  reported  in 
two  parts.  This  first  volume  briefly  discusses  the  models, 
but  concentrates  on  presenting  the  numerical  results  and 
analyses  of  the  acoustically  and  stochastically  modeled 
data.  Background  material,  references,  computer  program 
descriptions,  etc.  are  given  in  the  second  volume  in  the 
form  of  appendices.  The  contents  of  the  two  volumes  are 
summarized  below. 

Volume  I 

• Section  2 gives  a brief  description  of  the 

acoustic  models  used  in  the  study  to  generate 
signal  and  noise  data,  and  then  defines  the 
test  scenario  in  terms  of  the  ocean  environ- 
ment, target  geometry,  receiver  and  source 
characteristics,  etc.  Idealization  of  source 
and  receiver  properties  allows  the  investiga- 
tion to  concentrate  on  fluctuations  attributed 
to  a source  moving  through  multipath  inter- 
ference fields. 
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• Section  3 summarizes  the  basic  properties 

of  the  stochastic  process  models  to  be  tested: 
the  Gauss-Markov,  Jump,  and  Ehrenfest.  It 
also  describes  the  method  in  which  these 
models  are  used  to  simulate  fluctuations  in 
signal  and  noise. 

• Section  4 derives  and  compares  the  statisti- 
cal properties  of  acoustically  and  stochasti- 
cally modeled  transmission  los^  (TL)  fluctua- 
tions, which  themselves  determine  the  signal 
and  noise  fluctuation  characteristics.  The 
properties  investigated  are  limited  to  the 
first-order  statistics,  autocovariance,  and 
selected  level-crossing  properties,  all  of 
which  are  found  to  be  sensitive  to  the  defi- 
nition of  the  fluctuating  part  of  the  TL. 

The  density  function  and  spectrum  of  the 
acoustically  modeled  data  are  found  to  be 
consistent  with  the  theory  of  multipath 
propagation,  but  are  different  from  those 
asstimed  in  the  stochastic-process  models. 

• Section  5 considers  the  properties  of  acous- 
tlcally  modeled,  omnidirectional  ambient 
noise.  The  model  assiimes  that  surface  ships 
are  the  sources  of  noise,  and,  for  the  test 
scenario,  there  are  many  ships  contributing. 

As  expected  from  theory,  the  omnidirectional 
noise  is  usually  stable  (small  variance, 
large  decorrelation  times)  and  has  density 
function  consistent  with  that  for  many 
sources  of  multipath  fluctuations.  The 
exception  is  for  the  case  that  a single  ship 
passes  close  to  the  receiver  and  dominates 
the  noise  for  a long  time  period. 

• Section  6 examines  the  acoustic  and  stochas- 
tlc  ambient-noise  data  as  processed  by  the 
idealized  horizontal  towed  array.  Ten-hour 
time  series  of  acoustically  modeled  noise 
are  shown  to  be  sensitive  to  the  array's 
main-beam  width,  but  not  to  sldelobe  suppres- 
sion if  greater  than  30  dB.  As  expected  from 
theory,  the  noise  density  function  becomes 
broader  as  the  beamwldth  is  reduced  and  the 
variances  are  consistent  with  chi  square 
variables  whose  degrees  of  freedom  depend 
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on  the  nximber  of  ships  which  contribute  to 
the  fluctuations.  The  decorrelation  times 
were  large  and  not  very  sensitive  to  beam- 
width.  These  properties  of  the  acoustic 
data  were  found  to  be  reasonably  well- 
simulated  by  the  stochastic  processes  - 
since  the  marginal  densities  often  approached 
normal  near  the  median  and  the  autocovariance 
functions  decayed  smoothly  and  slowly.  As 
far  as  level-crossing  properties  are  con- 
cerned, the  Jump  process  showed  reasonable 
agreement  with  the  acoustic  data,  at  least 
for  lower  threshold  values. 

• Section  7 considers  time  series  of  signal- 
to-noise  ratio  (SNR)  as  generated  from  the 
acoustically  and  stochastically  modeled 
signal  and  noise  data  for  several  target 
tracks.  The  level-crossing  properties  are 
directly  associated  with  detection  histories, 
so  that  the  analysis  concentrated  on  them. 

For  the  best  stochastic-process  simulation 
of  signal  and  noise,  given  inputs  derived 
from  the  acoustic  data  themselves,  the 
detection-related  properties  were  found  to 
be  reasonably  well  simulated. 


Volume  II 


• Appendix  A provides  definitions  of  statisti- 
cal properties  for  stochastic  processes  and 
fluctuation  data  analysis. 

• Appendix  B gives  an  overview  of  sonar  per- 
formance prediction,  including  the  various 
types  of  acoustic  and  system  models  and 
levels  of  simulation. 

• Appendix  C describes  the  acoustic  signal  and 
noise  models  used  in  the  study,  and  discusses 
their  validity. 
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Appendix  D surveys  statistical  fluctuation 
models  based  on  an  acoustic  mechanism: 
varying  multipath  interference  caused  by 
source/receiver  motion.  The  various  statis- 
tical distributions  and  correlation  functions 
are  shown  to  be  closely  related. 

Appendix  E describes  and  lists  statistical 
properties  of  the  stochastic  fluctuation 
models  tested  in  the  study.  Although  not 
derived  from  properties  of  the  acoustic 
field,  these  models  are  in  common  use  in 
performance  modeling. 

Appendix  F discusses  the  approach  used  in 
the  study  to  simulate  a sonar  signal  pro- 
cessor and  detector.  Special  attention  is 
g paid  to  the  processor's  integration  times 

and  the  signal-excess  detector. 

Appendix  G provides  justification  for  cer- 
tain results  used  in  the  study  about  array 
response:  the  relative  Importance  of  trans- 
mission-loss and  beam-pattern-induced  fluc- 
tuations, conditions  under  which  beam  split- 
ting can  be  ignored,  and  the  effect  of  side- 
lobe  suppression  on  noise  fluctuations. 

Appendix  H documents  the  algorithms  used  to 
model  signal  and  noise,  the  numerical  pro- 
cedures for  generating  realizations  of  the 
random  processes,  and  the  statistical  ana- 
lysis package. 
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Conclusions  and  Recommendations 
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1.4.1  Conclusions 

On  the  basis  of  the  acousticallsr  modeled  data 
for  the  single  test  environment  and  limited  target/receiver 
parameters  of  this  study,  the  principal  conclusion  is: 
given  proper  input  data,  a stochastic-process  model  can 
produce  simulations  of  signal- to-noise  ratio  (SNR)  of 
sufficient  accxiracy  to  provide  useful  estimates  of  instan- 
taneous and  cumulative  detection  probabilities,  as  well 
as  expected  waiting  times  for  detection  and  loss  of  con- 
tact. In  addition,  the  random-process  models  considered 
here  were  found  to  be  capable  of  simulating  basic  proper- 
ties of  the  acoustic  time  series  of  signal  and  noise  with 
reasonable  accxiracy,  again  given  accurate  input  parameters. 
A careful  review  of  these  statements  is  in  order. 

First  of  all,  the  results  are  for  one  case, 
and  cannot  be  extended  without  some  further  analysis. 

The  implication  is  not  that  the  study  has  to  be  repeated 
in  detail  for  each  new  scenario,  but  rather  that  the  basic 
properties  of  other  acoustic  environments  and  SNR's  must 
be  somehow  estimated  (e.g.,  with  simplified  models)  and 
examined  for  amenability  to  stochastic-model  simulation. 

In  addition,  only  a few  statistical  properties  were  stud- 
ied; extrapolation  of  the  conclusions  to  Include  other 
detection-related  measures  cannot  be  made  with  confidence. 

Secondly,  the  stochastic-process  models  were 
given  the  best  possible  opportunity  to  agree  with  the 
acoustic  data;  the  detailed  trend  in  TL  and  the  calcu- 
lated statistical  properties  of  the  TL  and  noise  were 
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available  and  provided  the  inputs  to  the  models.  This 
is  not  the  usual  situation.  In  fact,  if  it  were  necessary 
to  run  acoustic  models  of  the  type  used  here  or  to  collect 
measurement  data  to  determine  the  stochastic  model  inputs 
for  every  new  case,  there  would  be  little  value  in  using 
the  stochastic  models  at  all:  the  "real"  time-series  data 
would  be  available. 

Finally,  although  the  basic  temporal  and  distri- 
butional properties  of  the  candidate  stochastic  processes 
were  generally  consistent  with  the  acoustically  modeled 
data,  they  were  not  in  good  agreement  - especially  in 
the  tails  of  the  distribution  and  the  form  of  the  auto- 
covariance function.  The  effect  may  not  be  apparent  in 
an  overall  review  of  detection  statistics,  for  it  is  in 
the  low  probability  events  that  such  details  become  impor- 
tant . 


In  spite  of  these  limitations,  the  conclusions 
encourage  the  use  of  the  random-process  approach.  In- 
cluded in  the  r'«coinmendatlons  given  below  are  suggestions 
for  removing  some  of  the  shortcomings  and  improving  the 
approach. 

1.4.2  Recommendations 

A substantial  amount  of  computer  machinery  has 
been  developed  or  applied  in  connection  with  this  effort: 
acoustic  noise  and  signal  simulation  models,  statistical 
analysis  and  detector  packages,  and  stochastic-process 
model  algorithms.  There  is  then  a capability  to  examine 
other  cases  and  other  problems  such  as: 
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- higher  acoustic  frequencies  and  short-range 
tactical  scenarios 

- detection  by  fields  of  sonobuoys 

- fixed  system  performance 

- correlation  for  multiple  arrays 

- more  complicated  detectors  and  tracking 
algorithms 

- higher  order  statistical  properties  of  the 
acoustic  parameters 

Given  the  analysis  routines,  models,  and  in- 
sight into  acoustic  fluctuation  mechanisms,  there  is 
also  the  opportunity  to  analyze,  interpret,  and  simulate 
measured  acoustic  fluctuation  and  performance  data. 

Although  "validation"  of  the  models  is  an  extremely 
ambitious  project,  careful  examination  of  a few  selected 
data  sets  from  well-controlled  experiments  would  enhance 
the  credibility  of  "model"  studies  such  as  this  one. 

There  is  a developing  foundation  of  theory  for 
multipath- induced  fluctuations  of  signal  and  noise  (Appen- 
dix D gives  a partial  survey).  It  can  be  applied  to 
the  present  problem  and  extended  in  several  ways; 

to  provide  input  parameters  for  the  stochas- 
tic-process approach  via  simplified  models 
or  analytical  techniques  for  predicting 
acoustic  fluctuation  properties  as  functions 
of  the  ocean  environment,  encounter  scenario, 
and  receiver  system. 
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- to  rank  order  the  scenario  parameters  accord- 
ing to  their  effect  on  fluctuations  and  detec- 
tion properties,  and  to  identify  conditions 
for  which  source-motion-induced  fluctuations 

% will  dominate  those  caused  by  the  dynamic 

ocean  medium 

- to  suggest  alternative  stochastic-process 
models  with,  e.g.,  more  realistic  marginal 
distributions  or  autocovariance  functions 

X 


As  an  example  for  the  last  point,  the  results  reported 
here  in  themselves  suggest  that  an  immediate  improvement 
to  the  stochastic  simulation  of  TL  fluctuations  for  the 
test  case  would  occur  if  the  Jump  process  were  to  use  a- 
chi-square  (or  non-central  chi-square)  marginal  distri- 
bution . * 


*k  stochastic  process,  derived  from  the  GM  process, 
with  this  marginal  distribution  has  in  fact  been  used 
in  TL  fluctuation  analysis  (see  Appendix  E). 


Section  2 

THE  TEST  CASE  AND  ACOUSTICALLY  MODELED  DATA 

Transmission  loss  and  noise  predictions  gener- 
ated from  state-of-the-art  acoustic  models  for  a special 
test  scenario  serve  as  the  control  data  for  this  study. 
They  provide  deterministic,  known  acoustic  environments 
which  can  be  viewed  as  "real  data"  for  the  test  of  se- 
lected stochastic  models.  They  also  offer  an  explanation 
for'  the  fluctuation  mechanisms  and  an  opportunity  to  im- 
prove the  stochastic  simulations. 

This  section  gives  brief  descriptions  of  the 
acoustic  models,  the  methods  used  to  generate  TL,  signal 
and  noise  data  for  the  specific  test-case  conditions, 
and  the  resulting  test  data  themselves. 

2 . 1 Acoustic.  Source,  and  Receiver  Models 

The  acoustically  modeled  data  are  generated  as 

follows . 


• Transmission  loss  (TL)  is  the  fundamental 
variable  and  is  predicted  as  a function 
of  range  with  a model  combining  Parabolic 
Equation  (PE)  TL  for  waterborne  paths  and 
FACT  TL  for  bottom  bounce.  The  environ- 
ment is  assumed  static,  so  that  the  pri- 
mary mechanism  for  fluctuations  is  the 
movement  of  a source  through  the  multi- 
path  interference  field.  For  the  environ- 
ment and  frequency  (25  Hz)  treated  here, 
the  model  is  fully  capable  of  predicting 
the  TL  to  the  detail  required  (the  range 
sampling  period  is  1/6  nm) . 
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• The  target  moves  at  constant  velocity  and 
radiates  a stable  tone  at  25  Hz. 

• The  receiver  is  an  idealized  towed  array, 
at  a fixed  location  and  depth  with  a sim- 
ple, deterministic  response  to  plane  waves 
(beam-pattern).  Note  that  the  use  of 
stable  source  and  receiver  characteristics 
allows  the  study  to  concentrate  on  the 
acoustic  fluctuations. 

• Signal  as  a function  of  time  is  calculated 
for  each  target  track  directly  from  the 
TL,  the  target  range  H(t)  and  bearing  (J)(t), 
the  source  level  SL,  and  the  array  response 
BP((j>)  according  to 

S(t)  - SL(t)  - TL(R(t))  + BP((|)(t)).  (2-1) 


• The  amhi ftnt.-noise  model  uses  the  same  equa- 
tion (2-1),  but  accounts  for  multiple 
^ sources.  It  uses  a hypothesized  set  of 

surface  ships,  moving  on  straight-line 
courses  across  the  basin  and  radiating 
noise  correlated  with  their  speeds.  Ship 
source  levels  and  courses  are  treated  as 
random  variables  subject  to  constraints 
% Imposed  by  somewhat  general  bounds  (e.g., 

speeds  are  15  knots  ±5  knots , courses  are 
within  15°  of  90°  or  180°T).  For  each 
point  in  time,  the  radiated  noise  from 
each  ship  propagates  to  the  receiving 
array.  The  resultant  contributions  are 
% added  incoherently  (random  phase),  under 

the  assumption  that  fluctuations  associated 
with  a coherent  sum  are  t3n?ically  averaged 
out  by  the  signal  processor's  time  inte- 
gration. Transmission  loss  is  derived 
from  the  model  described  above,  but  with 
t a source  depth  of  60  feet  and  corrected 

to  simulate  loss  from  a near-surface 
source  which  radiates  only  away  from  the 
surface.  The  simulated  noise  time  series 
is  sampled  every  two  minutes  for  the  dura- 
tion of  a replication  (e.g.,  10  hours). 

^ Each  replication  is  initialized  by  ran- 

dom selection  of  ships,  speeds,  courses 
and  source  levels . 
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More  detailed  information  on  the  approach  and  acoustic 
model  validity  can  be  found  in  Appendix  C.  Appendix  H 
describes  the  computer  routines  for  generating  signal 
and  noise. 

2 . 2 Test-Case  Scenario 


The  single  test  environment  is  in  deep  waters 
of  the  North  Pacific.  Not  only  does  it  exhibit  conditions 
of  operational  interest,  but  it  also  corresponds  to  a 
location  at  which  data  have  been  collected  and  the  TL 
model  evaluated.  The  ocean  environment  and  surface- 
ship  distributions  are  based  on  measvirements . 

The  test-case  conditions  are  summarized  below. 

2.2.1  Ocean  Environment 

The  sound  speed  profile  used  to  calculate  trans- 
mission loss  is  derived  from  measurement  data.  It  shows 
a thin  mixed-layer  at  the  surface,  a classic  thermocline, 
minimum  at  3000  feet , and  pressure  gradient  to  the  14000 
foot  bottom.  The  bottom  loss  is  also  based  on  measvire- 
ments, and  is  a constant  6 dB  per  bounce.  Independent 
of  grazing  angle.  Wind  speed  is  assumed  to  be  moderate 
and  thus  has  no  effect  on  the  calculations.  These  envlon- 
mental  conditions  are  assumed  to  obtain  throughout  the 
500-mile-radius  "basin."  The  discussions  on  model  val- 
idity of  Appendix  C support  this  assvimptlon. 
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2.2.2 


Transmission  Loss 


The  transmission  loss  calculations  to  500  nm 
are  made  for  a 1000  foot  receiver  and  for  sources  at  60 
and  300  feet.  The  model  assumes  narrowband  signals  at 
25  Hz.  The  TL  estimates  are  completely  deterministic 
and  computed  by  the  PE/FACT  model  mentioned  in  2.1. 
Transmission  loss  as  a function  of  range  is  sampled 
every  1/6  nm  and  used  for  both  signal  and  noise  fluctua- 
tions . 

2.2.3  Target 

The  target  is  assumed  to  have  a constant  source 
level  of  155  dB  (uPa)  at  25  Hz,  a constant  speed  of  5 
knots  and  a constant  depth  of  300  feet . Three  tracks 
have  been  used  for  this  study,  each  with  constant  course: 

Track  1:  Start  300  nm  due  north  of  the  receiver 
and  head  due  east 

Track  2:  Start  50  nm  due  north  and  travel 
due  east 

Track  3:  Start  300  nm  due  north  and  travel  due 
south 


Figure  2-1  shows  this  geometry. 
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2.2.4  Ship  Distribution 

The  shipping  distribution  used  in  the  calcula- 
tions of  the  noise  time  series  were  generated  by  the  DSBN 
model  mentioned  in  2.1  and  described  in  Appendix  H with 
the  following  parameters: 

• speed:  uniform  distribution  10-20  knots 

• bearing:  uniform  in  (75°,  105°)  and 

(255°,  285°) 

• average  inter-ship  arrival  time:  0.78  hrs 

• basin  size:  1000  miles  in  diameter,  array 

at  the  center 

• source  level:  161  + 60  log  (speed/15)  + 

20  log  [(160  + £)/152]  (dB 
re  lyPa,  Hz),  where  I is 
uniform  in  (-52,  +52). 

Once  determined,  each  ship  has  constant  course,  speed, 
and  source  level. 

This  model  is  consistent  with  the  traffic  pat- 
terns and  ship  types  observed  in  the  scenario  area.  The 
expected  ship  density  over  the  basin  is  one  ship  per 
4500  square  miles. 

2.2.5  Array 

The  array  model  described  in  Appendices  C and 
H was  used,  with  parameters  shown  in  Table  2-1.  The 
omnidirectional  response  is  also  calculated. 
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steering  Angle 

Beamwldth 
(half-power ) 

Sldelobe  Suppression 

90°  (Broadside) 

8° 

30  dB 

90° 

8° 

00 

90° 

4° 

30 

% 

90° 

4° 

00 

82°  (8°  east  of 
broadside) 

8° 

30 

74° 

8° 

30 

t 

. 65° 

8° 

30 

56° 

8° 

30 
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Array  Parameters 
Table  2-1 
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2.3 


The  Set  of  Acoustically  Modeled  Data 


% 


The  set  of  acoustically  modeled  data  generated 
for  this  study  with  the  models  of  2.1  and  the  scenario 
of  2.2  consists  of: 
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% 


% 


t 


% 


% 


t 


• Transmission  loss  versus  Range 

- sampled  every  1/6  nm  from  0 to  500  nm 

- for  receiver  depth  of  1000  feet  and 
source  depths  of  60  and  300  feet. 

- at  frequency  25  Hz 

• Signal  Level  versus  Time 

- sampled  every  2 minutes 

- for  target  tracks  1,  2,  and  3 

- for  source  depth  of  300  feet  and  re- 
ceiver depth  at  1000  feet 

- for  array  response  parameters  of  Table 
2-1 

- for  source  level  of  155  dB 

• Beam  Noise  versus  Time 

- 50  replications,  each  of  duration  10 

hours  and  one  replication  of  duration 
35  hours,  sampled  every  2 minutes 

- for  ship  parameters  described  in  2.2 

- for  an  omnidirectional  receiver  and 
for  the  array  response  parameters  of 
Table  2-1  (simultaneous) 

• Signal-to-Noise-Ratio  versus  Time 

- as  constructed  from  Signal  and  Noise 

time  series 

• Range-Averaged  TL 

- 2 nm  and  5 nm  intensity  averages  for 
above  cases 

- A+BlogR  fit  (LMS)  over  0-100  and  100-500  nm 
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Section  3 

THE  STOCHASTIC  MODELS  AND  THEIR  APPLICATION 


The  Gauss-Markov  (GM),  the  Gauss-Jump  (Jump), 
and  to  a limited  extent  the  Ehrenfest  random-process 
models  have  been  selected  for  testing  In  this  study,  pri- 
marily because  of  their  popularity  In  performance  and 
engagement  analyses . This  section  summarizes  their  ba- 
sic properties  as  well  as  the  manner  In  which  they  are 
used  In  this  study  to  simulate  fluctuations  In  signal, 
noise,  and  slgnal-to-nolse  ratio  (SNR). 

Additional  supporting  Information  can  be  found 
In  Volume  II.  Appendix  A gives  definitions  of  stochastic- 
process  properties.  The  rationale  for  the  choice  of  the 
selected  models  and  their  specific  properties  are  des- 
cribed In  Appendix  E.  The  computer  algorithms  for  gener- 
ating sample  functions  for  these  processes  and  the  proce- 
dures for  determining  Input  parameters  are  documented  In 
Appendix  H. 

3.1  Summary  of  Properties  of  Stochastic  Processes 

A stochastic  (or  random)  process  can  be  viewed 
as  a random  variable  which  depends  on  an  Index  (usually 
time).  The  three  models  used  In  this  study  are  mathe- 
matical In  nature,  with  little  basis  In  the  physics  of 
acoHStlc  propagation  and  noise  (although  each  Is  Inspired 
by  a physical  phenomenon).  Among  the  properties  common 
to  the  three  are; 
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• Each  depends  on  at  most  three  parameters. 

• Each  is  stationary  and  ergodlc  (for  the 
case  that  the  mean  is  zero). 

• Each  has  a one-dimensional  (marginal)  dis- 
tribution which  is  Gaussian  or  approxi- 
mately Gaussian. 


X 


X 


% 


t 


X 


t 


% 


Some  of  the  distinguishing  characteristics  of  interest  in 
this  study  are  listed  in  Table  3-1.  Noteworthy  are  the 
sample-path  structures  (see  Figure  4-5  for  examples  of 
continuous  and  Jump  paths)  and  the  exponentially  decaying 
autocovariance  function  which  yields  a smooth  spectrum  at 
low  frequencies,  with  no  dominant  components. 

3.2  Method  for  Simulating  Fluctuations 

The  random-process  models  are  used  to  simulate 
fluctuations  of  transmission  loss  with  range  and  of  sig- 
nal and  noise  with  time.  In  general,  let  Z(t)  denote 
one  of  these  quantities  and  let  X(t)  be  the  random- 
process  model  for  simulating  its  fluctuations.  The 
usual  approach  taicen  in  engagement  modeling 
the  variable  Z in  decibel  units,  to  predict 
minlstic"  part  Z(t),  and  then  to  simulate  Z 

Z(t)  - ?(t)  + X(t). 

If  Z is  viewed  as  a "trend,"  chosen  so  that 
fluctuation 

Y(t)  - Z(t)  - Z(t)  (3-4) 


is  to  treat 
its  "deter- 
with 

(3-3) 

the  residual 


t 
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GM Gauss- Jump Ehrenfest 


Parameters 

u,a,X 



n 

Index  Set 

tcro.Ti 

moBm 

k ■ 0,1,2, .. . 

Sample  Paths 

continuous , 
nowhere  dif- 
ferentiable 

piece-wise 
constant , 
discontinuous 
at  Jumps 

discrete 
(paths  con- 
structed by 
interpolation) 

Joint  Density 

2-dimen- 

sional 

Gaussian 

Not  Gaussian 

— 

Marginal 

Density 

Gaussian 

Gaussian 

(y,<j*) 

Binomial 

(normalized) 

Autoco- 

variance 

Function 

(n^)k 

Decorrelation 

Time 

1/X 

1/X 

-l/^n(S^) 

Power  Spec- 
trxim 

2a* 

2a* 

Level-Crossing  Properties  for  the  Gauss-Jump: 


• Expected  Waiting  Time  for  Jump  Process  to 
Exceed  L 

- (1/A)(1-Pj3)/Pj3  (3-1) 

• Probability  that  Jump  Process  Exceeds  L 
Sometime  In  Interval  [0,T] 

- l-(l-Pjj)exp(-Pj^T/X)  (3-2) 

where  P^  Is  the  probability  that  a Gaussian 
variable  ((i,<7*)  exceeds  L. 

Some  Properties  of  Stochastic  Processes 
Table  3-1 
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has  mean  zero,  then  the  random-process  model  should  have 

A 

mean  zero,  and  it  will  yield  a good  simulation  Z of  Z 
to  the  extent  that  the  properties  of  X agree  with  those 
of  Y. 


In  this  study,  the  trend  "ZCt)  of  the  acoustically 
modeled  series  Z(t)  is  found  by  averaging,  and  the  resi- 
dual fluctuation  series  Y(t)  is  analyzed  by  a statistical 
analysis  package  to  determine  its  sample  statistics  (dis- 
tribution function,  moments,  autocovariance  function, 
etc.).  These  statistics  of  Y are  in  turn  used  as  the  in- 
puts for  the  random-process  model  X. 

In  the  case  of  noise  time  series,  ? is  a con- 
stant, long-term  average.  For  transmission  loss  (TL(R)) 
versus  range,  the  trend  TL(R)  is  found  as  an  intensity 
average  over  a range  window  (2  nm  or  5 nm  or  the  entire 
interval  of  interest).  The  simulated  TL  is  then  given 
by 


TL(R)  « TL(R)  + X(R)  (3-5) 

The  associated  signal  time  series,  before  beamforming, 
is  derived  from  the  target's  range  R(t)  and  source 
level  SL(t)  as 


S(t)  - SL(t)  - [Ti:(R(t))  + X(R(t))].  (3-6) 

Finally,  note  that  the  stochastically  modeled 
data  are  derived  with  the  benefit  of  substantial  infor- 
mation about  the  acoustic  data  - viz. , the  trend  and  the 
statistical  properties  of  the  fluctuations.  In  the  usual 
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application,  only  the  trend  is  predicted  with  any  confi- 
dence (from  data  or  acoustic  models),  so  that  the  random- 
process  parameters  must  be  estimated  from  "experience." 

» Hence,  in  this  study,  the  stochastic  models  are  given 

the  best  possible  chance  to  generate  accurate  data. 
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RESULTS  FOR  ACOUSTICALLY  AM)  STOCHASTICALLY 
MODELED  TRANSMISSION  LOSS  AND  FLUCTUATIONS 


This  section  derives  and  compares  the  properties 
of  acoustically  and  stochastically  modeled  TL  and  its  fluc- 
tuations. Considerable  emphasis  is  placed  on  TL  since  in 
the  stochastic  data  it  is  an  essential  component  for  sig- 
nal, while  for  the  acoustic  data  it  drives  both  the  sig- 
nal and  noise  properties. 

Among  the  topics  covered  here,  one  of  the  most 
important  is  the  extraction  of  the  "trend"  or  "determinis-  i 

tic  part"  of  the  TL.  As  discussed  in  Section  3,  there 
are  a number  of  ways  in  which  this  can  be  done , and  the 
choice  impacts  heavily  on  the  properties  of  the  residual, 
which  is  then  the  fluctuation . 

The  first  part  of  the  section  presents  the  pro- 
perties of  the  TL  fluctuations,  and  relates  them  to  sig- 
nal variations.  In  the  second  part,  the  ability  of  the 
stochastic-process  models  to  simulate  these  properties 
is  investigated. 

4 . 1 Properties  of  Acoustically  Modeled  TL 

The  acoustic  TL  data  generated  for  this  study 
consist  of  a TL  range  series  generated  by  the  PE/FACT 
combination  model,  sampled  every  1/6  nm  from  0 to  500  nm, 
for  a frequency  of  25  Hz,  a 1000  foot  receiver,  and  two 
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source  depths.  The  analysis  here  concentrates  on  the 
data  for  the  deeper  (300  foot)  source  depth,  used  to  simu- 
late the  signal.  Figure  4-1  is  a plot  of  that  function. 

The  number  of  cases  (3  types  of  detrending  over 
several  range  intervals)  and  scope  of  the  study  preclude 
the  calculation  of  more  than  a few  properties  of  the  TL 
fluctuations.  The  ones  chosen  (see  Appendix  A)  are  those 
which  determine  stochastic-process  inputs  (mean,  variance, 
decorrelation  range),  those  which  can  be  easily  compared 
with  basic  stochastic-process  properties  (marginal  dis- 
tribution, autocovariance  function),  and  those  which  re- 
late to  signal  detectability  (TL  level-crossing  statis- 
tics) . 

4.1.1  Fluctuation  Properties  as  a Function  of  Smoothing 

The  usual  prediction  of  mean  TL  can  range  from 
a simple  curve  to  a detailed  average  over  short  ranges, 
and  the  resulting  fluctuation  properties  can  be  very 
sensitive  to  the  amount  of  smoothing.  To  investigate 
this,  an  average  of  the  TL(R),  called  TL(R),  was  con- 
structed using: 


(i) 

an  A+BlogR  fit 

(ii) 

a 5 nm  intensity 

average 

(iii) 

a 2 nm  intensity 

average 
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over  various  range  segments.  Then  the  "range  fluctua- 
tion" was  calculated  for  each  case  as 

ATL(R)  » TL(R)-fL(R). 

For  a source  moving  at  rate  R(t),  it  is  ATL(R(t))  which 
determines  temporal  multipath  fluctuations  and  it  is 
this  term  which  the  random-process  models  are  intended 
to  simulate.  Figure  4-2  shows  examples  of  ATL(R)  for 
the  three  types  of  TL  trends. 

4.1.2  First-Order  Statistics  of  ATL(R) 

The  first-order  statistical  properties  of  the 
various  versions  of  ATL(R)  have  been  calculated  and  are 
summarized  in  Table  4-1.  Notable  are  the  following: 
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• The  fluctuations  corresponding  to  the  de- 
tailed average  TL  (2  nm  or  5 nm  average) 
are  extremely  consistent  from  range-inter- 
val to  range-interval  beyond  about  50  or 
75  nm.  This  shows  an  absence  of  long- 
range  trends  (sometimes  used  as  an  indi- 
cator for  stationarity ) . This  is  not 

the  case  for  the  short-range  data,  which 
have  complicated  regions  of  bottom-bounce 
propagation;  nor  is  it  true  for  the  fluc- 
tuations about  the  A-^-BlogR  average,  which 
show  convergence  zones. 

• In  every  case,  the  positive  skewness  means 
the  ATL  distributions  are  skewed  to  the 
right:  the  long  tail  to  the  right  extends 
as  far  as  25  dB  and  corresponds  to  fades 
or  "drop-outs"  in  the  propagation,  while 
the  excursions  below  the  mean  are  limited 
to  about  5 dB. 

e Note  that  the  mean  is  not  necessarily  0 
when  TE  is  an  intensity  average,  and  in 
fact  is  positive. 
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Fluctuation  Statistics  for  ATL  • TL  - TL 


i 

1 


Now  in  studying  fluctuations  of  signal  inten- 
sity, it  is  -ATL  which  is  of  interest.  All  of  the  statis- 
tics of  Table  4-1  are  applicable,  with  the  following 
changes : 


The  sign  of  the  mean  and  skewness  are 
changed  indicating  skewness  to  the  left. 


The  percentiles  are  reversed  so  that  the 
long  tails  (up  to  25  dB)  are  to  the  left 
and  there  are  only  small  excursions  (5  dB) 
above  the  mean. 


The  resulting  properties  are  consistent  with  the  theory 
(Appendix  D)  of  intensity  as  a s\im  of  randomly-phased  multi- 
path  sinusoids.  In  fact,  the  various  candidate  intensity 
distributions  (chi-square,  noncentral  chi-square,  and  log- 
normal) were  transformed  into  dB's  and  tested  against  the 
sample  distributions  of  -ATL  over  the  100-500  nm  range 
interval.  The  results  showed  poor  agreement  for  the  log- 
normal, but  very  good  agreement  for  the  distributions 
predicted  by  theory  for  long-range,  waterborne,  multi- 
path  propagation  (chi-square  and  non-central  chi-square 
with  2 d. f . ) . 


4.1.3 


Temporal  Properties  of  TL  Fluctuations 


The  correspondence  of  range  fluctuations  in  TL 
with  temporal  fluctuations  of  signal  from  a moving  source 
has  been  emphasized  throughout  this  report.  In  Appendix 
D some  rationale  for  interpreting  the  spectral  properties 
of  the  intensity  in  terms  of  the  multipath  propagation 
structure  is  presented.  However,  the  scope  of  this  study 
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has  limited  the  investigation  to  one  of  the  temporal  pro- 
perties of  the  TL  (or  signal)  in  dB's,  as  used  in  the  sto- 
chastic models  for  signal  fluctuations.  In  particular, 
the  focus  is  on  the  autocovariance  function  form  and 
the  decorrelation  time  for  the  TL  fluctuations. 

Figure  4-3  shows  an  example  of  the  normalized 
autocovariance  functions  for  the  TL  fluctuations,  for  the 
three  TL  averages.  As  an  aid  to  identifying  the  periodi- 
cities, the  corresponding  Fourier  transforms  have  also 
been  calculated  and  samples  are  shown  in  Figure  4-4.  Im- 
portant features  derived  from  these  and  other  cases  in- 
clude: 
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% 
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• The  "decorrelation"  ranges  (to  the  1/6 
nm  resolution  available)  are  short  in 
every  case  for  fluctuations  about  the 
2 nm  and  5 nm  averages  - typically  1/3 
to  1/2  nm.  However,  for  the  A+BlogR 
case,  the  first  1/e  point  occurs  at  5 nm 
and  the  correlation  can  be  high  at  great 
ranges . 

• There  is  significant  structure  in  the 
autocorrelation  functions  beyond  the 
first  zero-crossing,  quite  different 
from  the  exponential  function  assumed 
in  the  stochastic  processes. 

• The  spectra  for  the  three  different 
averages  (^)  are  nearly  identical  above 
a cutoff  frequency  corresponding  to  the 
averaging  range.  For  the  100-500  nm 
range,  notice  important  contributions 

at  periods  from  about  1.5  to  4 nm  and 
near  30  nm,  with  the  last  eliminated  for 
the  case  of  detailed  averaging  (2  or  5 nm), 
as  expected.  For  the  shorter  range  inter- 
val , there  are  peaks  corresponding  to 
periods  in  the  neighborhood  of  0.5  to 
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TL  » A + BlogR 


Spectra  of  TL  Fluctuations 
Range  Interval  100-500  mn 
Figure  4-4 


2 nm,  plus  strong  low-frequency  components 
(periods  of  20'*’  nm) . These  two  scales  are 
consistent  with  the  theory  of  multipath 
interference  (Appendix  D)  for  waterborne 
paths  (small  vertical  angle  apertures)  and 
bottom  bounce  paths  (wider  aperture). 

4.1.4  Level-Crossing  Properties  of  TL 

Some  level-crossing  properties  of  the  total  TL 
(not  ATL)  are  summarized  in  Table  4-2  for  the  100-500  nm 
interval  and  the  300  foot  source  (per  Figure  4-1).  Note 
that  the  data  are  "inverted"  in  the  sense  that  for  sig- 
nal we  usually  look  for  properties  of  -TL  rather  than  +TL 
as  shown  here.  Also,  the  probabilities  in  the  table  are 
actually  sample  probabilities  estimated  from  the  2400 
data  points. 

If  the  three  levels  (99,  103,  107  dB)  are  viewed 
as  FOM's  and  if  the  detector  detects  whenever  the  TL  is 
less  than  the  FOM,  then  the  table  shows  that  for  a target 
with  any  range  in  100-500  nm  equally  likely,  there  are 
average  detection  probabilities  of  0.14,  0.44  and  0.74  for 
the  three  FOM’s.  The  expected  waiting  ranges  mean  that 
if  the  target  range  is  chosen  at  random  in  100-500  nm, 
and  the  target  is  opening,  then  the  waiting  time  to  cross 
the  level  L is  on  the  average  31,  1.8  or  0.3  nm,  respec- 
tively. The  small  standard  deviations  indicate  that  there 
is  not  much  variation  about  these  average  numbers. 

Table  4-3  displays  the  distribution  of  range 
Intervals  during  which  TL  is  always  below  (or  above)  a 


Number  of  Intervals  of  Length  k for  which  TL<L  (TL>L) 

Level-Crossing  Properties  of  Acoustically  Modeled  TL: 
Distribution  of  Interval  Lengths 
Over  Ranges  100-500  nm 

Table  4-3 
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level  L,  as  measured  over  the  400  nm  track  from  100-500 
nm.  The  interpretation  depends  on  the  way  in  which  inter- 
vals are  counted:  an  interval  of  length  k miles  is  one 
for  which  TL<L  (or  >L)  for  exactly  k miles  and  no  more. 

The  table  then  shows  the  number  of  intervals  within  a 
window  of  lengths.  For  example,  at  L » 99  dB,  there  are: 

2 intervals  with  lengths  between  0.16  and  0.2 
nm  for  which  TL<L 

9 intervals  with  lengths  between  0.5  and  1.0 
nm  for  which  TL<L 

- 6 intervals  with  lengths  between  0.5  and  1.0 

nm  for  which  TL>L. 


so  that  the  distribution  of  interval  lengths  is  constrained 
by  more  than  the  decorrelation  range. 


Some  observations  about  values  shown  in  the  table 

follow. 

• S3rmmetry  of  the  distributions  corresponding 
to  TL>103  and  TL<103  is  suggested.  This 
might  be  expected  since  P(TL>103)  is  near 
0.5.  However,  there  are  indeed  more  long 
intervals  with  TL>L  and  more  short  ones 
with  TL<L,  and  the  reason  is  that  the  TL 
fluctuations  are  skewed  - with  deep  fades 
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about  the  meau  (TL»L)  but  only  small 
excursions  of  TLi<L.  This  is  apparent 
in  the  waiting  times  of  Table  4-2  as  well. 


• Similarly,  differences  in  the  distributions 
of  interval  lengths  for  TL>99  and  TL<107, 
as  well  as  for  TL<99  and  TL>107,  are  attri- 
buted to  the  tendency  of  transmission  to- 
ward large  fades. 

• All  of  the  distributions  are  concentrated 
at  interval  lengths  between  0.5  and  2 nm. 
Although  a decorrelation  range  of  about  5 
nm  is  appropriate  for  the  total  TL  over  the 
100-500  nm  interval  (from  the  A+BlogR  fit), 
the  faster  fluctuations  seem  to  have  a sig- 
nificant influence  on  the  level-crossing 
properties . 

• Given  this  concentration  of  interval  lengths, 
the  principal  difference  between  one  distri- 
bution and  another  is  driven  by  the  number 
of  long  intervals.  Equation  (4-1)  shows 

how  a single  long  interval  displaces  many 
short  ones. 

These  features  will  become  important  when  the  stochastic 
data  are  compared  with  the  acoustic  data  in  the  next  sec- 
tion. 


4 . 2 Basic  Properties  of  Stochastically  Modeled  TL 

The  statistical  properties  of  the  stochastic  pro- 
cesses considered  in  this  study  were  summarized  in  Section  3. 
Since  each  process  is  ergodic,  after  a sufficient  amount  of 
time,  the  time-sampled  statistical  properties  will  approxi- 
mate the  ensemble  properties.  For  example,  if  the  process 
has  a normal  marginal  distribution,  then  a histogram  con- 
structed from  a sample  path  (realization)  will  approach 
that  distribution  as  the  duration  of  the  time  series  in- 
creases . 
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Even  though  most  of  the  properties  are  known, 
it  seems  worthwhile  to  summarize  the  results  for  40  repli- 
cations of  the  GM  and  Jump  processes  corresponding  to  the 
long-range  TL  fluctuations  of  4.1.  The  purpose  is  to 
compare  the  stochastic  and  acoustic  time  series,  and  to 
thus  illustrate  some  effects  of  finite  sampling  and  the 
differences  between  processes.  Along  with  the  TL  fluctua- 
tion properties,  the  level-crossing  statistics  of  TL  simu- 
lations (TL  + Random  Process)  will  be  displayed.  All 
of  this  is  used  to  judge  the  adequacy  of  the  stochastic- 
process  models  to  simulate  TL  and  its  fluctuations  as 
functions  of  average  TL. 

4.2.1  Sample  Functions  of  TL  Fluctuations 

As  an  aid  to  visualizing  the  differences  be- 
tween the  Jump  and  GM  processes.  Figure  4-5  shows  typical 
sample  functions  of  each.  They  correspond  to  the  para- 
meters (u,  <T,  t)  found  above  for  TL  fluctuations  over  the 
100-500  nm  interval,  and  for  the  2 nm  average  and  A+BlogR 
fit  for  TL.  For  the  A+BlogR  fit,  (u,  a,  t)  « (0.0, 

5.0,  4.7  nm),  and  the  two  sample  functions  clearly  exhi- 
bit the  special  character  of  their  parent  processes,  the 
step  versus  continuous  sample  paths.  The  number  of  jumps 
is  25  out  of  a total  of  1000  sample  points  for  the  Jump 
process  - consistent  with  the  value  of  t (one  jump  every 
4.7  nm  or  28  points).  To  the  more  detailed  ^ (2  nm), 
there  corresponds  a tighter  fluctuation  distribution, 

(u,  a,  t)  » (0.6,  2.2,  0.5  nm) , and  the  GM  and  Jiunp  pro- 
cesses begin  to  look  more  alike  (in  the  limit  of  t > 0, 
they  are  the  same). 


4-16 


Corresponding  to/  0 
TL  » 2 nm  Average  I 


5 *!!- 


Corresponding  to/  Q 
^ a 2 nm  Average  I 


Corresponding  to/  q 

TL  » A+BlogR  I 

(-5 


».uiijni. 1 ■■■■«■  ■««fii 


Corresponding  to/  Q 
TL  ■ A+BlogR  I 


50  nm 


100  nm 


150  nm 


Figure  4-5  Sinulated  TL  Fluctuaticns 


Comparison  of  the  sample  paths  of  Figure  4-5 
with  those  of  Figure  4-2  suggests  that  the  GM  process 
better  mimics  the  character  of  the  acoustic  simulations, 
but  that  for  t small  (detailed  TL)  the  Jump  process  may 
also  be  quite  satisfactory. 

4.2.2  Sample  Moments  and  Distribution  Functions 

Sample  mean,  standard  deviation,  skewness,  and 
kurtosis  were  calculated  for  the  simulated  TL  fluctuations 
described  in  the  preceding  paragraph  and  are  displayed  in 
Table  4-4.  The  "long-term"  and  ensemble  values  of  these 
statistics  are  known,  but  it  is  instructive  to  notice 
that  the  subinterval  moments  make  the  (stationary)  GM 
and  Jump  processes  look  even  more  non-stationary  than 
the  acoustically-modeled  fluctuation  process.  This  is, 
of  course,  a vagary  of  sampling:  the  observed  values 
are  within  expected  limits  corresponding  to  the  number 
of  independent  samples  making  up  the  statistics.  (The 
decorrelation  times  shown  in  the  same  figure  are  used 
to  estimate  the  "effective"  number  of  "independent" 
samples).  Ensemble  statistics  were  also  calculated  for 
40  replications  of  the  stochastic  processes;  they  look 
much  the  same  as  the  time  statistics  and  are  not  displayed. 

The  sample  distribution  functions  were  calcu- 
lated for  40  replications  of  the  processes.  As  expected, 
when  the  number  of  Independent  samples  is  properly  esti- 
mated (using  decorrelation  and  zero-crossing  times),  the 
sample  functions  are  accepted  at  the  0.05  or  higher  sig- 
nificance levels  as  normal  (l.e.,  log-normal  intensities) 
with  a Lilliefors  version  of  the  Kolmogorov  test.  This 
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- .3 

- .3 

.7 

4.7 

- .0 

3.0 

300-500 

-3.1 

4.4 

0.0 

O.l 

-1.1 

3.9 

- .3 

0.3 

100-500 

-1.1 

4.7 

- .1 

-.3 

-0.3 

4.4 

- .4 

3.0 

^Decorrelaticn  Range  and  Range 

3.3 

to  First  Zero-Crossing  of  Auto- 

4 

covariance  (nm) 

3.9 

19 

Effective  Nuxiner  of  ” Indepen- 

dent"  Points  Over  100-500  ran 
* (out  of  2400  points) 

100 

130 

Ol  Corresponding 
to  2 nm  Average  TL 

(u»0.6;  0-2.2;  t-0.5  nm) 


JUMP  Corresponding 
to  2 nm  Average  TL 

(u^.6;  0-2.2;  tM).5  nm) 


Statistics  for  a TVpical  TL  Sinulation 
Over  Range  Interval  100-500  m 


Table  4-4 
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is  of  course  one  major  discrepancy  between  the  acoustic 
and  stochastic  simulations  - since  the  acoustic  TL  fluc- 
tuation distributions  were  rejected  as  normal  for  the 
detailed  TL.  We  note  also  that  the  stochastic  process 
simulations  often  passed  the  Kolmogorov  test  at  the  0.01 
level  for  the  high  d.f.  chi-square  intensity  distribu- 
tion. The  reason  is  that  any  Gaussian-like  haystack  den- 

\ 

sity  function  looks  log-normal  near  the  mean. 

4.2.3  Correlation  Properties 

Normalized  autocovariance  functions  for  the 
four  sample  functions  described  above  are  shown  in  Figure 
4-6.  Again,  note  that  the  GM  and  Jump  samples  are  typi- 
cal of  the  many  which  were  generated  for  this  study. 

The  plots  show  the  expected  exponential  decay.  The  ob- 
served decorrelation  ranges  are  near  those  used  as  . 
input  to  the  random  process  generator  (see  Table  4-3). 

The  sample  covariance  functions  become  "noisy"  for  lags 
beyond  about  10%  of  the  total  range  of  observation  (400 
nm) . 


The  difference  in  correlation  properties  be- 
tween the  stochastic  and  acoustic  simulations  is  apparent 
when  the  curves  of  Figure  4-6  are  compared  with  those 
of  Figure  4-3.  However,  a better  appreciation  can  be 
gained  when  the  corresponding  Fourier  transforms  ( "spec- 
tra") of  Figures  4-7  and  4-4  are  compared.  The  spectra 
for  the  stochastic  processes  are  as  expected  for  an  ex- 
ponential autocovariance  function:  they  have  the  proper 
shape  (from  Table  3-1): 

S(a»)  - 0)-  2irf,  X - l/x. 
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Autocovariance  Functions  for  Simulated  TL  Fluctuations 
♦ (Over  Range  Interval  100-500  nm) 

Figure  4-6 
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These  three  range-periods  are  about  30,  15,  and  10  nm  for 
the  TL  » A+BlogR  case  (t  = 5 nm)  and  about  3,  1.5,  and 
1 nm  for  the  2-mile-average  case  (t  = 0.5  nm).  The  spec- 
tra for  the  acoustically  modeled  TL  fluctuations  (Figure 
4-4)  are  not  at  all  similar  to  these,  and  the  same  obser- 
vation holds  for  simulations  over  other  range  intervals. 

This  second  discrepancy  between  properties  of 
the  two  types  of  models  can  be  summarized  as  follows. 


• Whereas  the  stochastic  simulations  display 
exponentially  decaying  autocovariance  func- 
^ tions  and  their  monotone,  smooth  spectra, 

the  acoustically  modeled  data  show  signifi- 
cant structure  in  the  autocovariance  func- 
tions (resembling  a damped  cosine)  and 
spectra  with  "power”  concentrated  at  pre- 
ferred fluctuation  periods  (e.g.,  1 to 
I 5 nm ) . 


4.2.4 


Level -Crossing  Properties  for  TL(R) 


The  level-crossing  statistics  of  the  stochastic 
TL  simulations  were  obtained  from  samples  of  simulated 
'fL(  R)  generated  as 
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'^(R)  - TL(R)  + A^L(R) 

over  100<R<500  mn  with  TL  either  (A+BlogR)  or  a 2-mile 

A 

% average  and  ATL(R)  the  appropriate  GM  or  Jump  process 

(i.e.,  the  same  cases  focussed  on  throughout  this  section). 

A 

For  40  replications  of  the  400  nm  range-series  of  TL,  the 
statistics  were  calculated  for  comparison  with  those  of 
♦ the  acoustic  TL(R)  (displayed  in  Table  4-2).  These  level- 

crossing statistics  are  summarized  in  Tables  4-5  and  4-6 

with  the  20^^,  50^^,  and  80^^  percentile  values  of  P(TL<L), 

A 

P(TL>L),  and  expected  waiting  "time.s  " In  the  same  tables 
^ are  displayed  the  corresponding  acou.  .ically  modeled  sta- 

tistics (from  Table  4-2). 


• The  results  for  TL  » (A+BlogR)  show  a wide 
range  of  values  from  the  20th  to  80th  percen- 
tile, and  the  best  agreement  is  for  the  J\imp 
process . 

• For  the  2-mile  average  the  range  of  values 
is  much  smaller  and  the  Jump  process  again 
shows  best  agreement  with  the  acoustic  data. 

In  fact,  except  for  the  waiting  time  for 
TL<99  dB,  the  median  values*  are  very  close 

to  the  statistics  for  the  acoustic  data. 


t Table  4-7  shows  distributions  of  interval  lengths 

for  which  TL  is  less  than  L.  These  are  based  on  "typical" 
replications,  and  in  fact  the  general  properties  of  these 
statistics  were  quite  consistent  from  one  sample  function 
% to  another.  The  corresponding  data  for  the  acoustically 

modeled  TL  (from  Table  4-3)  are  given  in  the  same  table, 
and  a comparison  shows: 

*Note  that  the  median  values  are  in  good  agreement  with 
I the  analytic  level-crossing  solutions  for  the  Jump  pro- 

cess with  a constant  threshold  (very  nearly  the  case 
for  TE  - A+BlogR),  as  given  by  formulas  (3-1)  and  (3-2) 
of  Table  3-1. 
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Simulation 


Interval  Length  k(nm) 

0.16-0.2  i 0.2-0. 5 i 0. 5-1.0  ! 1-2 ! 2-4  I 4-10 i 10-16  ; >16 


L - 99  dB 


Log  Fit  Jump 
2-mlle  Avg  + Jump 
Log  Fit  -f  Of 
2-oille  Avg  + GM 


107  dB  Log  Fit  + Jump 

2-nlle  Avg  -f  Jump 
Log  Fit  + Q! 
2-mlle  Avg  + QI 

Acouatlc 


Number  of  Intervals  of  Length  k with  TL<L 

Level-Crossing  Properties  of  Acoustically  Modeled  TL: 
Distribution  of  Interval  Lengths 
Over  Ranges  100-500  nm 

Table  4-7 
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• The  GM  and  Jump  processes  are  readily  dis- 
tinguishable by  their  level-crossing  re- 
sults for  short  interval  lengths.  For 
intervals  greater  than  the  decorrelation 
range,  they  show  similar  behavior, 

• The  GM  process,  for  either  smoothed  TL, 

consistently  overestimates  the  number  of  ‘ 

short  intervals. 

• The  Jump  process  added  to  TL  =•  A+BlogR 
generally  overestimate  the  number  of  long 
intervals . 

• The  best  match  over  all  cases  ^ given  by 
the  sum  of  the  2-raile  average  TL  and  the 
Jump  process.  The  agreement  is  not  very 
good  - and  this  seems  to  be  attributable 
to  the  lack  of  concentration  of  interval 

lengths  between  0.5  and  2 nm.  In  light  of  ^ 

equation  (4-1)  and  the  discussion  of  4.1.4, 
the  random-process  models  with  their  expo- 
nential autocovariance  functions  and  S3nn- 
metric  distribution  functions  yield  sample 
paths  with  too  many  short  intervals  and  i 

too  few  medium-length  intervals. 

i 

In  conclusion,  the  level-crossing  properties  of 
the  acoustically-modeled  TL  over  a long  range  interval 
(100-500  nm)  are  simulated  moderately  well  by  the  Jump 
process  added  to  a 2-mile  intensity-averaged  TL.  No 
other  combination  of  TL  and  ATL  is  a very  close  second. 

Another  test,  over  shorter  range  intervals  but  for  signal- 
to-noise  ratios,  is  given  in  Section  7. 


4 . 3 Summary 

Acoustically  modeled  TL  as  a function  of  range 
(R)  is  treated  as  a sum  of  two  components:  the  trend, 
fL(R),  and  the  fluctuation,  ATL(R)  - TL(R)  - tE(R). 
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Statistical  properties  of  ATL(R)  for  various  range  inter- 
vals within  0 to  500  nm  have  been  calculated  for  three 
types  of  smoothing:  2-mile  intensity  average,  5-mile 
average,  A+BlogR  fit. 

The  GM  and  Jump  processes  were  used  to  simulate 
ATL(R),  with  input  parameters  (mean,  variance,  decorrela- 
tion time)  determined  directly  from  the  acoustic  data. 

A 

Sample  functions  of  ATL(R)  were  also  added  to  TL  to  yield 

A A 

simulations  of  TL:  TL  * TL  + ATL. 

Highlights  of  the  comparison  of  the  acoustically 
and  stochastically  modeled  TL  are; 

(i)  The  sample  ma.rgi’^al  density  functions  for 
the  acoustic-intensity  fluctuation  data 
match  a non-central  or  central  chi-square 
density  better  than  the  log-normal  of  the 
stochastic  process. 

(ii)  The  sample  autocovariance  functions  for 

the  acoustic  fluctuation  data  showed  struc- 
ture different  from  the  exponential  form 
of  the  stochastic  models.  Fourier  trans- 
forms illustrate  that  fluctuation  power 
occurs  at  preferred  frequencies,  in  con- 
trast to  the  stochastic-process  models' 
smooth,  monotone  spectra. 

(iii)  Selected  level-crossing  properties  of  acous- 
tically modeled  TL  itself  were  best  fit  with 
a simulation  using  a 2-mile  average  ^ and 


the  Jump  process.  However,  the  agreement 
was  only  moderately  good,  with  differences 
possibly  attributable  to  (i)  and  (ii). 

For  signal  modeling,  the  detection-related  properties  of 
TL  are  its  level-crossing  statistics.  Hence,  for  the  case 
at  hand,  the  Jump-process  fluctuation  (given  the  proper 
inputs)  added  to  the  2-mile  average  TL  can  yield  a satis- 
factory simulation.  The  impact  of  inaccuracies  in  the 
simulation  on  performance  predictions  is  discussed  further 
in  Section  7. 
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Section  5 


ACOUSTICALLY  MODELED  OMNIDIRECTIONAL  AMBIENT  NOISE  DATA 


Fifty  ten-hour  and  one  thirty-five-hour  repli- 
cations of  the  acoustic  noise  model  for  the  conditions 
given  in  Section  2 were  run  for  the  study  of  beam-noise 
modeling  and  detection.  One  of  the  byproducts  is  a set 
of  omnidirectional  noise  time  series,  and  it  is  profitable 
to  examine  these,  both  for  comparison  with  the  beam  noise 
data  given  below  and  for  verification  that  sidelobe  levels 
are  as  modeled  in  Appendix  G.  No  stochastic  simulations 
of  omni  noise  were  tested. 

On  average  there  were  155  surface  ships  within 
the  500  nm  radius  which  defines  the  ocean  basin.  Because 
of  this  large  number  of  contributing  sources,  the  noise 
is  expected  to  be  quite  stable  over  time,  and  the  statistics 
bear  this  out.  Figures  5-1  and  5-2  show  time-series  re- 
cords for  two  typical  replications.  In  the  first  (called 
"Case  1"),  a ship  passes  close  to  the  receiver  and  per- 
turbs the  very  stable  field  for  several  hours.  The  corre- 
sponding autocovariance  function  is  nearly  linear  with  a 
decorrelation  time  (t)  of  about  one  hour.  The  spectrum 
shows  only  that  the  noise  is  dominated  by  a low-frequency 
trend  (period  of  at  least  several  hours).  The  mean  is 
83.7  dB  and  a « 2.9  dB.  The  second  case  shows  a number  of 
perturbations  by  nearby  or  noisy  ships,  but  no  large 
excursions.  The  result  is  a mean  of  83  dB,  but  a much 


Time  Series,  Autocovariance,  Spectrum  of  Typical  Omni  Noise  Sample  (#1 ) 

Figure  5-1 
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smaller  standard  deviation  (1.2  dB).  The  autocovariance 
function  and  spectrum  of  Figure  5-2  have  a little  more 
character,  reflecting  the  semi -periodic  occurrence  of 
slightly  higher  levels  every  hour  or  so.  Although  the 
decorrelation  time  is  about  20  minutes  (as  compared  to  an 
hour  in  Case  1),  the  first  zero-crossing  of  the  autocovari- 
ance function  is  beyond  two  hours  here,  as  opposed  to  about 
1.5  hours  in  the  first  case. 


The  large  values  of  the  decorrelation  times  indi- 
cate that  moments  or  distribution  functions  obtained  by  time 
sampling  over  ten  hours  will  be  poor  estimates,  there  may 
be  as  few  as  5 or  10  "independent"  samples.  As  a partial 
solution,  we  assume  ergodicity  and  calculate  moments  at 
several  fixed  time  points  over  the  ensemble  of  50  repli- 
cations. The  results  are  extremely  consistent  from  time 
point  to  time  point; 


Percentile 

10 

50 

90 

mean 

83.3  dB 

83.7 

84 

standard  deviation 

1.5  dB 

2.1 

2.8 

skewness 

.96 

1.7 

2.9 

kurtosls 

0.6 

2.1 

4.2 

The  implication  is  that  the  omni  noise  distribution  has 
small  variance  and  positive  kurtosls.  The  skewness  to 
the  right  is  obviously  the  result  of  the  occasional  appear- 
ance of  high  noise  levels  from  nearby  or  loud  ships.  In 
spite  of  the  likelihood  of  under samp ling,  the  sample  den- 
sity functions  for  Cases  i and  2 above  are  shown  in  Figure 
5-3,  together  with  the  densities  for  log  transforms  of 
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exponential  and  chi-square  (15  d.f.)  variables.  The  latter 
is  a "best  fit"  for  Case  2.  While  Case  1 is  "bimodal" 

(as  expected),  Case  2 could  pass  as  log-normal  or  chi- 
square,  consistent  with  the  theories  presented  in  Appen- 
dix D. 

Also  note  the  variation  from  replication  to 
replication  of  the  decorrelation  times  of  the  omni  noise 
times  series: 


Decorrelation  Time  (mins) 
First  Zero-Crossing  (mins) 


Obviously,  questions  about  long-term  trends  in  noise  will 
require  much  longer  sample  time  series  than  are  considered 
in  this  study. 


There  is  one  long  replication  of  noise  (35  hours), 
and  it  is  encouraging  to  note  that  its  covariance,  spec- 
triim  and  statistics  (Figure  5-4,  Table  5-1)  are  consistent 
with  those  of  the  10  hour  series.  There  is  indeed  little 
correlation  beyond  the  first  zero-crossing  time  and  the 
long  term  statistics  show  the  c.haracteristic  positive 
skewness. 


In  summary,  omnidirectional  ambient  noise  for 
a case  in  which  many  ships  contribute  is  quite  stable 
(a  over  ten  hours  is  t3rpically  less  than  3 dB)  with  first 
order  properties  consistent  from  sample  to  sample  (the 
mean  ranged  from  83  to  84  dB,  the  skewness  from  1 to  3). 
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• Decorrelation  Time  * 56  minutes 

% 

• Time  to  First  Zero-Crossing  of  Autocovariance 
■ 156  minutes 


Statistics  for  35-Hour  Omni  Noise  Time  Series 

Table  5-1 
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The  decorrelation  times  are  large  and,  in  view  of  the 
small  variance,  not  too  important.  It  is  observed  that 
the  principal  uncertainty  in  predicting  such  noise  pro- 
perties is  in  estimating  the  time  and  duration  of  nearby 
ship  "interference."  Simple  analytical  models  which 
might  eliminate  that  problem  are  available. 


Section  6 


% 


♦ 


♦ 


> 


* 


RESULTS  FOR  ACOUSTICALLY  AND  STOCHASTICALLY 
MODELED  BEAM-NOISE  FLUCTUATIONS 

This  section  derives  and  compares  the  properties 
of  acoustically  and  stochastically  modeled  beam  noise,  i.e., 
ambient  noise  as  observed  by  the  generic  towed  array  system  I 

described  in  Subsection  2.2.  j 

I 

j 

Before  the  results  are  presented,  certain  sim-  s 

I 

plifications  in  the  analysis  should  be  noted.  The  array  I 

response  is  ideal,  represented  by  a time-invariant  beam  I 

pattern  with  a simple  main  lobe  and  uniformly  suppressed  | 

sidelobes  (see  2.2).  As  discussed  in  Appendix  G,  the 
choice  of  array  steering  angles  near  broadside  makes  multi- 
path  beam  splitting  unimportant , while  the  sidelobe  sup- 
pression (30  dB  or  more)  guarantees  that  noise  sources  off 
the  main  beam  contribute  little  to  the  fluctuations. 

Finally,  recall  that  in  the  study  of  transmis- 
sion loss,  a predictable  trend  was  removed  and  the  resi- 
dual became  the  fluctuating  component.  In  the  case  of 
noise,  however,  only  the  long-term  (e.g.,  24  hours)  mean 
is  assumed  to  be  predictable,  so  that  any  departure  from 
the  mean  is  treated  as  a fluctuation. 

6. 1 Properties  of  Acoustically  Modeled  Beam  Noise 

The  acoustic  beam-noise  data  were  generated 
with  the  DSBN  model  described  in  Appendices  C and  H,  and 


6-1 


I 


t 


¥ 


♦ 


I 


♦ 


¥ 


used  the  TL,  shipping  and  array-response  inputs  identified 
in  Subsection  2.2. 

Fifty  realizations  of  DSBN  beam  noise,  each  of 
duration  ten  hours,  were  generated  for  this  study.  The 
model  was  run  with  independently  selected  ship  parameters 
each  time,  and  a few  sample  cross-correlations  suggest 
that  the  50  series  can  be  treated  as  independent  samples. 
The  time  series  were  themselves  sampled  at  two  minute  in- 
tervals so  that  each  realization  consists  of  300  time 
points.  For  each  of  the  ten-hour  DSBN  series,  outputs 
were  obtained  for  eight  array  beams  (Table  2-1)  as  well 
as  the  omnidirectional  noise.  There  is  also  one  long  (35 
hour)  time  series  of  noise,  so  that  long-term  trends  can 
be  investigated  for  at  least  one  case. 

6.1.1  Typical  Beam-Noise  Time  Series 

Begin  by  looking  at  two  typical  noise  time 
series,  the  same  two  discussed  in  Section  5 as  Case  1 
and  Case  2.  Figures  6-1  and  6-2  show  the  time  series  for 
near-broadside  beams  with  8°  beamwidth  and  30  dB  side- 
lobe  suppression.  Although  the  time  scales  resemble  those 
for  omnidirectional  noise  (Figures  5-1  and  5-2),  the  beam- 
noise  data  show  much  larger  excursions.  The  reason  is 
that  there  are  on  average  155  ship  sources  contributing  to 
the  orani  noise,  but  only  seven  ships  for  the  beam  noise  (8° 
beam  and  ambiguous  image  beam). 

The  complicated  behavior  of  beam  noise  can 
be  appreciated  when  the  simultaneous  outputs  of  three 
different  beams  (Case  2)  are  compared  as  in  Figure  6-3. 
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Time  Series,  Autocovariance,  and  Spectrum  for  Beam  Noise  (Case  it2) 
90”  Beam  (0°  T),  8°  Beamwidth,  30  dB  Sidelobe  Suppression 


The  plots  show  the  effects  of  ships  moving  this  way  and 
that  across  the  array.  It  is  also  instructive  to  compare 
with  the  Omni  noise  of  Figure  5-2. 
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In  each  of  the  sample  series  notice  that  the 
fluctuations  are  rather  slow  - i.e.,  the  periods  are  long. 
Because  the  shipping  lanes  are  nearly  perpendicular  to 
the  beam  axis,  typical  radial  ship  speeds  are  on  the  order 
of  0 to  5 knots.  Since  the  important  TL  fluctuation  periods 
are  on  the  order  of  1 nm  or  more,  the  fastest  fluctuations 
should  have  periods  of  about  ten  or  fifteen  minutes.  Like- 
wise, the  time  for  a ship  to  cross  the  beam  is  on  the  or- 
der of  one  hour  at  100  nm  and  five  hours  at  500  nm.  These 
figures  are  consistent  with  the  pictures.  Also  notice 
that  as  the  steering  angle  decreases,  the  low  amplitude 
fluctuations  become  more  rapid,  showing  that  the  TL  fluc- 
tuations become  more  important  as  noise  sources  cross  the 
beam  with  greater  radial  speeds. 

6.1.2  Sidelobe  Contributions 

The  effect  of  sidelobe  suppression  is  typified 
by  Figure  6-4  and  the  upper  plot  of  Figure  6-3.  Here 
noise  from  the  broadside  beam  with  4^  beamwldth  and  30  dB 
sidelobe  suppression  is  compared  with  that  for  infinite 
suppression.  From  the  discussions  of  Appendix  G,  the 
greatest  sensitivity  is  expected  at  narrow  beamwidths. 

Even  in  this  case,  the  differences  are  apparent  only  at 
the  lowest  levels  of  noise.  Statistics  of  the  fluctua- 
tions, including  the  autocovariance,  for  fifty  replica- 
tions showed  that  there  is  no  appreciable  difference 
between  the  30  dB  and  Infinite-suppression  time  series 
of  beam  noise. 
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Time  Series,  Autocovariance  and  Spectrum  for  Beam  Noise  (#2) 
90^^  Beam.  4*^  Beamwldth,  30  dB  Sldelobe  Suppression 
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6.1.3  Moments  and  Distributions 

Table  6-1  shows  percentiles  for  the  sample 
mean,  standard  deviation,  skewness,  and  kurtosls  derived 
from  the  fifty  ten-hour  noise  time  series  and  various 
beams  and  beamwldths.  Since  there  is  probably  an  under- 
sampling  problem  (there  may  be  as  few  as  ten  uncorrelated 
samples  in  a ten  hour  series),  the  ensemble  statistics 
over  the  fifty  replications  have  also  been  calculated 
at  15  time  points;  the  results  are  displayed  in  Table 
6-2. 


The  first  thing  to  notice  is  that  the  ensem- 
ble and  time  estimates  of  mean  and  variance  are  nearly 

identical,  supporting  the  assumption  of  ergodicity.  Other  | 

significant  properties  of  these  data  include: 


1 


♦ 


♦ 


(i)  The  mean  levels  decrease  by  about  4 dB 
as  the  beamwldth  decreases  from  8^  to 
4°.  At  the  same  time,  the  variance  in- 
creases by  at  least  1 dB  and  the  skew- 
ness tends  to  decrease. 

(ii)  The  levels  increase  by  about  1.5  dB  from 
the  90°  to  the  65°  beam  (8°  beamwldth), 
although  the  variance  and  skewness  seem 
to  be  invariant. 

(iii)  Refer  to  the  omnl  noise  data  (Section  5), 
and  note  that  the  8°  beamwldth  levels 
are  down  about  17  dB  from  the  mean  omnl 


Beam 

Width 


(8°) 


(8®) 


(9°) 


(8°) 


(8°) 


Sldelohe 

Suppres- 

sion 


30  dB 


30  dB 


30  dB 


30  dB 


30  dB 


Standard  | 

Deviation  I Skewness 


65.8  dB  3.4  dB 
67.0  4.1 

68.5  5.0 


66.0  dB  2.7  dB 


67.0 

67.8 


67.5  dB  2.8  dB 

68.4  3.8 

69.5  4.5 


65.6  dB 

66.9 

68.5 


61.0  dB 
62.6 
64.3 


66.2  dB  4.3 


*The  three  numbers  in  each  box  are  the  20^^,  50^^,  and  80^*^ 
percentiles  for  50  ten-hour  replications. 


Time-Sampled  Statistics  of  Beam  Noise 


Table  6-1 
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O 00 

30  dB 

1 

65.7  dB* 

66.7 

67.8 

1 

3.8  dB 

4.3 

4.7 

65.3  dB 

2.4 

82° 

30  dB 

66.0 

2.8 

(8°) 

68.3 

4.5 

67.8  dB 

3.4 

65° 

30  dB 

68.4 

3.9 

(9°) 

69.0 

4.4 

(8°) 


30  dB 


65.5  dB 

3.9  dB 

66.3 

4.5 

67.5 

4.8 

61.6  dB 

5.4  dB 

62.9 

5.7 

63.6 

6.2 

*20^^,  50^^,  and  80^^  percentiles  over  50  ten-hour  repli- 


cations . 


Ensemble  Statistics  of  Beam  Noise 
(15  Time  Points  Over  50  Replications) 
Table  6-2 
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levels,  while  the  standard  deviation 
for  the  8°  beanrwidth  is  twice  that  for 
the  Omni  case.  The  skewness  of  the 
Omni  data  is  significant  and  positive 
(1-3)  while  that  for  the  beam  data  is 
typically  near  zero  or  negative. 

(iv)  The  moments  for  the  35-hour  sample  are 
consistent  with  the  other  statistics. 

The  only  point  which  requires  further  discussion  is  the 
change  in  levels  from  omni  to  8°-beam  noise  to  4°-beam 
noise.  Since  the  expected  number  of  ships  contributing 
in  each  case  decreases  from  155  to  8 to  4,  the  mean 
intensity  is  expected  to  decrease  by  about  13  dB  and 
3 dB  Instead  of  the  17  dB  and  4 dB  observed.  This  dis- 
crepancy is  ascribed  to  the  fact  that  the  means  derived 
from  the  time  series  are  averages  of  dB-variables . which 
have  different  skewness  according  to  the  number  of  ships 
contributing.  The  4 dB  and  1 dB  differences  result  from 
the  displacement  of  the  mean  of  the  log  from  the  log  of 
the  mean.  Similarly,  the  increase  by  1.5  dB  in  mean 
level  from  the  broadside  (90°)  beam  to  the  65°  beam  is 
attributed  to  an  Increase  in  the  expected  number  of  ships 
(from  eight  to  nine)  and  the  change  in  the  distribution 
of  ships  present. 

Sample  distribution  functions  for  besun  noise 
were  constructed  from  the  time-series  data  but  were  typi- 
cally undersampled.  As  the  small  positive  and  negative 
skewness  of  the  beam  noise  would  suggest,  the  sample 


(intensity)  distributions  were  most  often  accepted  by 
the  Kolmogorov  test  as  log-normal  or  chi-square  (2  d.f. 
and  higher)  at  the  5%  level. 

6.1.4  Correlation  Properties 

The  autocovariance  functions  for  beam  noise 
shown  in  Figures  6-1,  6-2  and  6-4  are  representative 
examples  of  the  50  replications.  They  are  adequately 
fitted  by  decaying  exponentials. 

Table  6-3  gives  the  percentiles  of  decorrela- 
tion times  for  the  various  beams  as  derived  from  the  .50 
replications.  The  most  concrete  conclusion  to  be  made 
is  that  both  beam  noise  and  omni  noise  have  decorrelation 
times  of  0.5  to  1.0  hour  and  zero-crossings  at  1.5  to 
2.5  hours.  The  fluctuation  time  scales  for  ships  on  the 
beam  (periods  of  at  least  10  or  15  minutes  driven  by  TL) 
and  the  long  waiting  periods  for  ships  to  arrive  on  a 
beam  (.8  hours)  are  consistent  with  these  figures. 

The  spectra  of  Figures  6-1,  6-2,  and  6-4  rein- 
force the  argument  that  the  long-term  trends  (periods  of 
one  hour  or  more)  can  dominate  the  temporal  properties 
of  both  beam  noise  and  omni  noise.  On  the  other  hand, 
when  a few  noise  sources  have  significant  radial  velocity 
components,  the  faster  (TL-induced)  fluctuations  also 
become  important  - as  noted  in  the  case  of  the  65*^  beam. 
Appendix  G gives  rules  of  thumb  for  determining  the  rela- 
tive importance  of  beam-pattern  and  TL  fluctuations  for 
a single  moving  source. 
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Beam-Noise  Decorrelation  Times 
(50  Replications  of  Ten-Hour  Duration) 
Table  6-3 
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6,1.5  Level-Crossing  Properties 

Although  noise  replications  were  generated 
for  nine  array  beams,  the  level-crossing  properties  are 
calculated  only  for  the  broadside  beam  with  8°  beamwidth 
and  30  dB  sidelobe  suppression.  The  statistics  are 
summarized  in  Table  6-4  (the  first  number  in  each  tri- 
plet) for  threshold  levels  equal  to  the  mean,  u ■ 67.1  dB, 
and  to  u i cr.  They  were  derived  from  many  ten-hour  reali- 
zations of  the  DSBN  model  and  represent  average  values. 

The  spreads  about  them  were  wide  in  all  categories.  The 
nearly  symmetric  and  Gaussian  nature  of  the  data  is  re- 
flected in  the  values  P(N>ii)  and  P(N<L).  The  cumulative 
probabilities  and  waiting  times  are  consistent  with  the 
average  decorrelation  time  of  36  minutes. 

Finally,  for  the  distribution  of  interval 
lengths  with  N>L  or  N<L,  the  case  N<L  is  focussed  upon 
here.  Two  samples  from  the  fifty  replications  of  beam 
noise  were  selected  as  representative  of  the  diversity  in 
the  distribution  of  lengths.  They  might  be  viewed  as  20th 
and  80th  percentile  curves,  since  they  were  chosen  with 
that  in  mind,  although  no  objective  criteria  were  applied. 
Table  6-5  shows  counts  of  interval  lengths  for  three  thres- 
hold levels  ( L * n,  u-*’a,  u-a)  and  the  two  samples  under 
the  labels  Acoustic  Simulation  #1,  #2.  The  values  are 
representative  of  data  which  consistently  show  up  to  two 
short  intervals  (<10  minutes),  one  or  two  intermediate- 
length  intervals  (10-60  minutes),  and,  according  to  L, 
zero  to  three  long  intervals  (60-200  minutes)  in  a 600 
minute  sample.  A single  long  interval  at  the  lowest  thres- 
hold level  occurs  when  the  beam  sees  only  a few  distant 
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Acoustic 
Acoustic  ^2 


GM  #1 
GM  #2 
Jump  #1 
Jump  #2 


Acoustic  #1 
Acoustic  f?2 


GU  #1 
GM 

Jump  #1 
Jump  #2 


Acoustic  ^1 
Acoustic  ^2 


GM  #1 
GM  #2 
Jump  #1 
Jump  M2 


Number  of  Intervals  of  Length  k with  N<L 

Level-Crossing  Properties  for 
Representative  Samples  of  Acoustic  and  Stochastic 
Beam-Noise  Time  Series  of  Ten-Hour  Duration 

Table  6-5 
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ships  or  ships  between  convergence  zones.  The  probability 
of  such  a situation  enduring  for  an  hour  or  more  in  a ten 
hour  period  is  high  for  the  8°  beamwidth  and  the  given 
ship  distribution. 

6 . 2 Properties  of  Stochastically  Modeled  Beam  Noise 

The  moments,  sample  distributions,  autocovari- 
ance, and  spectra  for  the  GM  and  Jump  processes  have  all 
been  discussed  in  Subsection  4.2.  Although  the  GM  pro- 
cess looks  more  like  the  beam-noise  series,  the  Jump 
matches  the  above-mentioned  properties  just  as  well. 

In  fact,  the  acoustic  simulations  show  one-dimensional 
distributions  which  are  nearly  Gaussian  and  autocovari- 
ance functions  which  are  approximately  exponential. 

As  far  as  the  level-crossing  characteristics 
are  concerned.  Table  6-4  compares  svimmary  data  for  the 
random  processes  with  those  for  the  acoustic  model  dis- 
cussed above  (90°  beam,  8°  beamwidth,  etc.).  Except  for 
the  cumulative  probabilities  and  the  wait  for  NMj  at 
L - p + a,  the  two  processes  yield  estimates  which  (on 
average)  are  in  good  agreement  with  the  acoustic  noise 
data.  There  is  again  an  undersampling  problem  Inherent 
in  these  numbers  since  there  are  in  each  ten-hour  series 
an  average  of  only  about  20  "uncorrelated"  time  samples. 
Nonetheless,  the  aggregated  results  suggest  that  the  ran- 
dom process  models  are  reasonable  candidates  for  simu- 
lating beam-noise  data. 

*For  example,  the  chance  of  having  eight  or  fewer  ships 
on  the  main  beam  (and  its  image)  and  all  beyond  200  nm 
is  more  than  0.1. 
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For  a final  measure,  compare  the  distribution  of 
interval  lengths  (with  N<L)  for  sample  paths  of  the  ran- 
dom-process model  wich  the  acoustic  data  in  Table  6-5. 
Just  as  for  the  acoustic  simulation  distributions,  the 
values  labled  "GM  #1,"  ”#2,"  and  "Jump  #1,"  "#2"  were 
derived  from  representative,  but  diverse  samples  paths 
of  the  GM  and  Jump  processes.  They  can  also  be  compared 
with  the  entries  of  Table  4-7  for  the  A+BlogR  average 
TL,  which  are  calculated  on  a different  scale  but  again 
representative.  Comparison  of  values  for  these  samples 
suggests  the  conclusions  that ; 

• As  noted  in  the  simulation  of  TL,  the  GM 
process  tends  to  predict  too  many  short 
intervals  and  too  few  long  ones. 

• The  Jump  process  is  in  reasonably  good 
agreement  with  the  acoustic  data  for  all 
values  of  L. 

6 . 3 Summary 

For  the  test  cases  considered  here,  the  Jump 
and  GM  random-process  models  can  adequately  simulate  the 
acoustic  beam-noise  data,  provided  that  accurate  inputs 
(p,a,T)  are  used.  This  conclusion  is  similar  to  that  for 
the  transmission-loss  data;  however,  the  estimation  of 
input  parameters  may  be  even  more  difficult. 

As  noted  in  6.1.3,  the  mean,  variance,  and 
shape  (e.g.,  skewness)  of  the  beam-noise  distributions 
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are  quite  sensitive  to  beam  width;  they  and  the  decorre- 
lation time  do  not  necessarily  depend  on  the  shipping 
field  in  a simple  way.  Calculating  these  parameters  with 
a large-scale  computer  simulation  such  as  DSBN  eliminates 
the  value  of  the  random-process  model  - since  "real"  noise 
time  series  and  their  properties  are  then  available. 

As  alternatives,  the  chi-square  and  shot  models 
mentioned  in  Appendix  D provide  analytic  techniques  which 
can  potentially  yield  adequate  estimates  of  the  shape  of 
the  distribution  and  autocovariance  functions.  In  fact, 
the  dependence  of  u,a  and  skewness  on  beamwldth  observed 
in  the  test  cases  can  be  explained  with  a chi-square  mo- 
del when  the  number  of  degrees  of  freedom  are  proportional 
to  the  number  of  contributing  ship  sources. 

Finally,  the  beam-noise  time-series  plots  show 
that  there  can  be  significant  excursions  in  levels  caused 
by  the  arrival  of  a loud  or  nearby  ship  on  the  beam  or 
by  the  absence  of  ships  on  the  beam.  The  narrower  the 
array's  main  lobe,  the  more  important  are  these  occurrences. 
The  details  of  system  performance  depend  on  these  "beam 
free  times"  and  durations  of  higher  or  lower  noise.  The 
shot  model  provides  a method  for  simulating  such  ship 
movements  across  array  beams. 
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Section  7 


RESULTS  FOR  ACOUSTICALLY  AND  STOCHASTICALLY 
MODELED  SIGNAL-TO-NOISE  RATIO 
AND  DETECTION  PROPERTIES 

To  this  point,  transmission  loss  and  noise  have 
been  considered  separately,  and  it  was  concluded  that  for 
the  test  cases  the  stochastic  models  can  simulate  statis- 
tical properties  of  the  acoustic  data  reasonably  well, 
provided  that  accurate  input  parameters  (2-mile  average 
TL;  (u,<7,t)  for  ATL  and  noise)  are  used.  In  this  section, 
the  TL  and  beam  noise  are  combined  to  form  the  signal-to- 
nolse  ratio  (SNR),  the  variable  which  determines  whether 
or  not  detections  are  made. 

For  a stable  source  with  level  SL,  the  SNR  as 
a function  of  time  is  given  approximately  by* 

SNR(t)  « SL  - TL(R(t))  + BP(<Kt))  - N(t) 


where 


R(t)  is  the  target  range  at  time  t, 

(^(t)  is  the  target  bearing  at  time  t, 

BP  is  the  array  intensity  response  (in  dB's), 

N(t)  is  the  beam-noise  level  at  time  t. 

Thus,  for  the  SNR  detector  assximed  in  this  study,  the  level 
crossing  properties  of 


*For  ship  speeds  much  smaller  than  sound  speed 
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t 
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» 


t 


BP((0(t))  - TL(R(t))  - N(t) 

determine  detection  properties.*  The  relative  importance 
of  the  variables  BP,  TL,  and  N as  fluctuation  components 
has  been  discussed  in  Appendix  G and  will  be  noted  for  the 
test  cases  below.  Clearly,  if  R(t)  is  nearly  constant 
(target  distance  changes  slowly)  then  BP  will  drive  the 
signal  fluctuations.  On  the  other  hand,  if  the  target 
moves  on  a radial  path  from  the  receiver,  then  the  range 
change  will  be  rapid,  BP  constant,  and  TL  fluctuations 
will  dominate  for  signal. 

This  section  begins  with  examples  of  SNR  time 
series  to  illustrate  the  Interaction  of  the  signal  and 
noise  components.  Comparisons  of  the  detection-related 
statistics  for  the  acoustic  and  stochastic  simulations 
are  then  given  for  the  test  cases. 


' 7.1  Examples  of  Acoustically  Modeled  Time  Series 

^ of  SNR 

Consider  the  examples  of  acoustically  modeled 
time  series  of  slgnal-to-nolse  ratio  (SNR)  at  the  array 
% beamformer  output  shown  in  Figures  7-1  and  7-2.  Without 

the  benefit  of  prior  information  about  the  TL  and  noise 
properties,  one  might  note  some  trends  in  the  data,  but 
rely  heavily  on  "randomness"  to  explain  the  variations 
t in  SNR.  However,  the  plots  use  a single  replication  of 

noise  (Case  #1  of  Sections  5 and  6)  and  the  detailed 


♦The  iSnA  detector  has  a constant  recognition  differential 
R^ , and  the  signal  excess  (SE)  is  then  SE(t)  * SNR(t)-R^. 
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TL  of  Figure  4-1.  The  only  difference  then  among  plots 
are  the  target  tracks: 

• The  "tangential"  target  in  Figure  7-1  is 
opening  at  5 knots  on  a course  perpendicu- 
lar to  the  broadside  beam  (Track  1 of  Figure 
2-1) 

• The  "radial"  target  in  Figure  7-2  is  closing 
the  receiver  at  5 knots  on  a track  along  the 
axis  of  the  broadside  beam  (Track  3 of  Figure 


% 


X 


I 

* 

I 

I 

I 

I 


Radial 

Target 


_L 

Receiver 


■*  Tangential 
Target 


Concentrate  first  on  the  tangential  target  (Figure 
7-1).  Since  the  radial  range  rate  is  only  about  0.7  knots, 
the  signal  level  is  nearly  constant  as  the  target  crosses 
the  beam  (3  dB-down  points  for  array  response  are  shown 
as  dashed  lines)  and  the  SNR  variations  are  those  of  the 
beam  noise.  In  fact,  the  details  of  the  plot  of  beam  noise 
(Figure  6-1)  can  be  found  (inverted)  in  Figure  7-1. 

The  radial-target  case  is  more  complex  since  it 
is  Influenced  by  both  the  TL  and  beam-noise  fluctuations. 
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The  data  of  Figure  7-2  show  the  slow  but  pronounced  noise 
variations  observed  above  in  the  tangential-target  case. 
However,  the  noise  modulates  the  more  rapid  and  regular 
signal  fluctuations.  For  the  short  range  case  (100-50  nm), 
the  signal  has  a pronounced  period  of  about  10  minutes, 
corresponding  to  the  TL  range  period  of  about  1 nm  identi- 
fied in  Section  4.  When  the  target  is  more  distant  (300- 
350  nm),  the  range  period  is  on  the  order  of  1.5  nm,  so 
that  signal  fluctuation  periods  are  closer  to  15  or  20  minutes. 

7 . 2 Comparison  of  Models  for  the  Tangential-Target  Case 

In  the  tangential-target  case,  the  properties  of 
beam  noise  determine  the  detection  properties.  The  In- 
fluence of  the  beam  pattern  on  the  signal  is  of  course  im- 
portant, but  once  the  target  is  on  a beam  (e.g.,  within 
the  3-dB-down  bearings),  it  is  minimal.  Detection  proper- 
ties at  various  detector  threshold  levels  can  thus  be  de- 
termined for  both  acoustic  and  stochastic  models  from  the 
noise  level-crossing  properties  discussed  in  Section  6 
for  N<L*. 


The  comparison  of  fifty  replications  of  the  acous- 
tic SNR  with  many  sample  paths  for  the  stochastic  processes 
leads  to  the  conclusions  for  this  case  that : 

e For  thresholds  within  4 dB  of  the  mean  SNR, 

the  Jump-process  simulation  of  SNR  (actually, 
of  beam  noise)  was  Judged  to  be  in  reasonably 
good  agreement  with  the  acoustic  data. 


♦¥o  convert,  a noise  level  of  67  dB  corresponds  to  SNR  ■ 
-15  dB. 
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0 The  GM-process  results  were  not  as  good, 

but  still  satisfactory.  The  Ehrenfest  model 
was  also  used  in  this  test,  and  showed  re- 
sults about  the  same  as  the  GM. 

Note  well  that  the  input  parameters  for  the  random-process 
models  (noise  ctatistics  and  mean  TL)  were  determined 
directly  from  the  acoustic  data.  -Without  the  advantage 
of  accurate  Inputs,  the  stochastic  simulations  can  be  ex- 
pected to  fall. 

Consider  next  a few  of  the  comparisons  on  which 
these  conclusions  are  based.  Since  the  detection  parame- 
ters are  the  quantities  ultimately  used  in  performance 
analyses,  the  test  statistics  are  presented  in  those  terms. 
Table  7-1  compares: 

• detection  probabilities  (Pjj)  for  detectors 

requiring  the  SNR  to  exceed  the  threshold 
continuously  for  2,  10,  and  30  minutes 

e cumulative  detection  probabilities  (CPq) 

for  60-mlnute  opportunities 

e expected  waiting  time  for  the  2-mlnute  de- 
tector to  gain  contact 

e expected  waiting  time  for  the  2-mlnute  de- 
tector to  lose  contact 

e all  for  mean  signal  excess (^)  of  -4,  0, 
and  +4  dB. 
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The  table  entries  are  derived  from  many  replications  of 
the  acoustic  noise  model  and  the  Gauss-Jiunp  model  (using 
inputs  derived  from  the  acoustic  data).  Except  for  one 
waiting-time  entry,  the  Jvimp  process  shows  good  agreement 
with  the  acoustic  data.  The  small  discrepancies  in  the 
other  values  can  be  explained  by  the  Jump-model's  use  of 
a symmetric  (normal)  density  for  noise  and  the  observed 
skewness  (to  the  right)  of  the  density  for  the  acoustic 

noise  data. 

7.3  Comparison  of  Models  for  the  Radial-Target  Case 

The  radial-target  case  is  more  complex  and  chal- 
lenges the  stochastic  models  since  it  is  strongly  influenced 
by  TL  fluctuations.  Compare  the  acoustically  modeled  SNR 
of  Figure  7-2  with  the  stochastic  data  in  Figure  7-3.  The 
latter  is  based  on  the  difference  of  two  Jump-process  sam- 
ple paths,  one  for  TL  and  one  for  beam  noise.  In  particu- 
lar, the  TL  for  signal  is  simulated  as  the  sum  of  the  2- 
mlle  average  ^ and  a Jump  process  whose  Inputs  (u,  a,  t) 
are  derived  from  the  acoustic  TL  fluctuation  data  over 
2-150  nm  for  the  nearby  target  and  over  100-500  nm  for  the 
distant  one.  Likewise  the  noise  is  simulated  as  a second 
Jump  process  with  (y,  a,  t)  taken  from  the  acoustic  beam- 
noise  data.  Although  the  plots  of  Figures  7-2  and  7-3 
show  similarities,  the  nulls  and  periods  are  Indeed  differ- 
ent in  character.  But  the  important  properties  are  mea- 
sured by  the  detection  statistics. 

As  in  the  tangential-target  case,  replications 
of  the  acoustic  SNR  were  generated  and  compared  to  sample 
paths  for  various  combinations  of  the  random-process  models. 
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• For  the  test  cases , the  random-process  mo- 
dels for  simulating  TL  and  beam  noise  for 
both  near  and  distant  radial  target  can  be 
combined  to  yield  adequate  simulations  of 
signal-to-noise  ratio,  provided  that  accurate 
input  parameters  are  used.  Especially  impor- 
tant is  the  smoothed  TL  (a  5-mile  or  finer 
average  is  indicated)  and  mean  beam-noise 
level , 

• Of  the  stochastic  model  combinations  consi- 
dered, the  Jump  model  for  both  TL  and  beam- 
noise  fluctuations  is  Judged  to  give  the 
best  overall  agreement ' with  the  acoustic 
results. 

Some  examples  of  the  comparisons  leading  to  these  conclu- 
sions follow. 

Tables  7-2  and  7-3  show  the  same  statistical  pro- 
perties for  the  radial  target  tracks  as  were  displayed  in 
Table  7-1  for  the  tangential  target.  The  values  of  recog- 
nition differential  in  these  cases  (R^j  » -6,  -12,  and  -18 
dB)  were  chosen  to  give  a reasonable  spread  of  conditions 
of  mean  signal  excess.  In  each  table,  the  results  of  the 
acoustic  simulation  are  compared  to  averages  derived  from 
the  stochastic  model  described  above  (J\imp  for  both  TL  and 
noise  fluctuations),  judged  to  be  the  best  of  those  consid- 
ered. 
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arlson  of  Acoustically  and  Stochastically  Modeled  Detection  Properties 
5-knot  Radial  Target  on  Broadside  Beam  of  Width  8°,  at  Range  300-250  ni 


For  the  nearer  target  (Table  7-2),  the  principal 
discrepancies  are  the  overly  optimistic  predictions  by  the 
stochastic  model  of  P^'s  and  Waiting  Times.  Just  as  in 
the  tangential-target  case,  these  may  well  result  from  the 
difference  in  shapes  of  the  marginal  densities  for  the 
acoustic  data  and  the  (Gaussian)  random-process  densities. 

For  the  distant  target  (Table  7-3),  the  discrepan- 
cies are  small,  and  not  clearly  due  to  any  single  property 
of  the  stochastic  model.  The  lack  of  complicated  bottom- 
bounce  propagation  is  one  reason  for  better  behavior  of 
the  data  for  this  case. 

7.4  Summary 

To  compare  acoustic  and  stochastic  simulations  of 
SNR  and  detection  properties,  this  study  has  concentrated 
on  two  special  target  tracks.  In  one  case,  the  target 
maintains  a nearly  constant  range  to  the  array,  and  crosses 
the  broadside  beam  perpendicularly.  Since  the  TL  changes 
very  little,  the  beam  response  and  noise  properties  deter- 
mine the  SNR  fluctuations.  In  the  second  case,  the  tar- 
get moves  radially  toward  the  array,  along  the  broadside 
beam  axis.  The  TL  fluctuations  are  then  important,  and 
combine  with  the  beam-noise  variability  to  yield  a more 
complicated  picture.  These  two  cases  are  bounding  in  the 
sense  that  for  any  other  target  track  the  TL,  beam-pattern, 
and  noise  fluctuations  will  occur  in  combination,  somewhere 
in  between  the  two  (Appendix  G offers  some  n^les  of  thumb 
to  estimate  the  relative  importance  of  TL  as^  array-response 
fluctuations) . 


For  both  the  tangential  and  radial  target  cases, 
the  best  agreement  with  acoustically  modeled  detection  pro- 
babilities was  obtained  for  the  stochastic  simulation  con- 
sisting of 


• Jump  model  for  beam  noise 

• Jump  model  for  TL  fluctuations  and  a 2-mlle 
averaged  TL 

The  agreement  was  Judged  to  be  reasonably  good  over  a range 
of  thresholds,  but  depends  critically  on  the  use  of  accurate 
Input  parameters: 

e Smooth  TL  and  mean  noise  level 
e TL  and  beam-noise  fluctuation  statistics 

(m,  <7,  T). 


Distribution  List 


Address 


1.  Defense  Advanced  Research  Projects 
Agency 

Washington,  DC  20301 

(Tactical  Technolo^  Office) 
(Technical  Library) 


2.  Assistant  Secretary  of  the  Navy  (R,E&S) 
ATTN:  Dr.  D.  Hyde 
The  Pentagon 
Washington,  DC  20350 

3.  Defense  Documentation  C^ter 

Cameron  Station  \ 

Alexandria,  VA  22314  \ 

t 

4.  Chief  of  Naval  Operations  \ 

Department  of  the  Navy  y 
Washington,  DC  20350  / \ 

(OP-987)  / \ 

(OP-955)  ' 

5.  Chief  of  Naval  Kpterlal 
Department  of  the  Navy 
Washington,  IXT  20360 

(ASW;^ 

6.  Offlcf^f  Naval  Research 
80^.  Quincy  Street 
A/lIngton,  VA  22217 

(Code  431) 

(Code  222) 

(Code  102-0S) 

(Code  486) 


7.  Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20360 

(Code  06H2) 

8.  Naval  Research  Laboratory 
Washington,  DC  20375 

(Code  8100) 


9.  Naval  Electronics  Systems  Conmand 
Department  of  the  Navy 
Washington,  DC  20360 
(PME-1 24-62) 

(NAVELEX  320) 


Distribution  List 


Address  Copies 

10.  Center  for  Naval  Analysis 
1401  Wilson  Boulevard 

Arlington.  VA  22209  1 

11.  Naval  Ocean  Systems  Center 
San  Diego.  CA  92152 

ATTN:  H.  Schenk  1 

12.  Naval  Underwater  Systems  Center 
New  London  Laboratory 

New  London,  CT  06329 

ATTN:  Mr.  R.  Masse  1 

13.  Bolt,  Beranek  and  Newman,  Inc. 

1701  North  Fort  Myer  Drive 
Arlington,  VA  22209 

ATTN:  Or.  M.  Moll  1 

14.  Operations  Research,  Inc. 

1400  Spring  Street 
Silver  Spring,  MO  20910 

ATTN;  Dr.  E.  Moses  1 

15.  Daniel  H.  Wagner  Associates 
Station  Square  One 

Paoll,  PA  19301 

ATTN:  Dr.  B.  McCabe  1 

16.  Systems  Control,  Inc. 

1801  Page  Mill  Road 
Palo  Alto,  CA  94304 

ATTN:  Dr.  J.  Anton  1 

17.  Institute  for  Acoustic  Research 
615  Southwest  Second  Avenue 
Miami,  FL  33130 

Mr.  J.6.  Clark  (ATTN)  1 

18.  Naval  Ocean  Research  andDevelopment  Activity 
Code  320 

Bey  St.  Louis,  MS  39520  1 


